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FOREWORD  J 


This  Is  the  final  report  covering  work  performed 
for  the  Naval  Surface  Weapons  Center  under  contract 
N60921-75-C-0007 , The  purpose  of  the  contract  was  to 
analyse  simplified  propagation  codes  and  recommend 
Improved  models  for  characterising  the  propagation  of  high 
energy  CW  laser  beams.  The  work  was  performed  under  the 
direction  of  Mr.  Larry  Jobson  of  NSWC. 

This  report  documents  the  work  performed  during 
the  period  January-Merch  1975.  K also  provides  a review 
of  the  earlier  work  (July  1974-January  1975).  The  earlier 
work  was  described  in  detail  in  interim  reports  SAI-74-587- 
WA  (October  1974)  and  SAI -74-622-WA  (January  1975). 

The  authors  wish  to  express  their  appreciation  to 
Dr.  Tom  Tuer  of  the  SAI  Ann  Arbor  office  for  his  develop- 
ment of  the  multi-line  absorption  coefficients.  We  also 
thank  Mrs.  T.  Peckham  for  implementing  the  numerous  code 
modi fications , 
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Section  X 
INTRODUCTION 


In  support  of  the  Systems  Analysis  Team  of  the 
Naval  Suffice  Weapons  Center,  Science  Applications,  Inc. 
(SAI),  performed  a series  of  tasks  designed  to  develop 
a bettor  understanding  of  the  modeling  a 1 ternat i ves 
available  to  laser  systems  analysts  and  to  assist  in  the 
continual  improvement  of  the  HOI  Engagement  Code  (NOtEC). 
Emphasis  was  placed  on  the  resolution  of  propagation 
Issues  of  interest  to  systems  analysts.  The  tasks 
involved  quantitative  comparison  of  simplified  HEL 
propagation  codes  used  by  DOD  analysts,  clari f icat i on  of 
the  beam  quality  problem,  development  of  beam  distortion 
and  displacement  models,  development  of  a simplified 
optical  train  model,  and  expansion  and  improvement  of 
the  AFWt  COMBO  code. 

The  work  was  divided  into  six  tasks.  The  first 
four  tasks  are  reviewed  in  this  report  (Section  4)*  but 
are  covereo  more  extensively  in  interim  reports  (References 
1 and  ?): 

First  I n t erim  Report  - Qul y-Septemh or  1974  , 

OTT5T 5 B }~WA7 'October  1974  ( 1' e f e"r o nee  lT” 

Task  1.  Quantitative  Comparison  of 
Simplified  Propagation  Coder' 

Task  3,  Beam  Quality  Modeling. 

Second  Interim  Report  - September  1P75- 
vfanu" a" r v T 9T5' , S A f-T4 - 6T2 - W A , January  T5? $ 

TOTe  cence 

Task  2.  Beam  Shape  and  Displacement 
Due  to  Thermal  tensing. 

Task  4.  Effect  of  Truncation  and 
Obscuration  on  the  Far -Field  Beam  Profile. 


Task  $ required  the  development  of  a simplified 
optical  train  model.  This  work  is  described  in  Section  2 
of  this  report.  Task  6 involved  modification  and  extension 
of  the  AFWI  CONDO  code  and  included  \he  incorporation  of 
models  developed  in  the  earlier  tasks.  The  new  version 
of  the  code  was  named  SAtCOM.  Section  3 describes  the 
modifications  and  provides  a user’s  guide  to  SAiCON. 
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Section  2 

A SIMPLIFIED  OPTICAL  TRAIN  MODEL 


When  performing  laser  application  studies,  the 
system  analyst  all  too  frequently  assumes  that  the 
characteristics  of  the  beam  leaving  the  transmitting 
optic  are  the  same  as  the  characteristics  of  the  beam 
leaving  the  laser  device'  without  regard  for  any  changes 
induced  by  the  optical  elements  and/or  components  required 
to  get  the  beam  from  the  laser  device  to  the  transmitting 
optic,  i.e.,  the  optical  train.  It  is  k;own,  however, 
that  significant  changes  to  the  characteristics  of  the 
beam  do  occur  as  it  propagates  through  the  optical  train 
(see  for  example  References  3 and  4).  The  finite  absorp- 
tivity of  the  high  power  mirrors,  clipping  and  blockage 
of  a portion  of  the  high  power  beam,  and  diffraction 
effects  all  reduce  the  available  power  in  the  beam  at  the 
transmitting  optic.  In  addition  the  optical  quality  of 
the  beam  is  affected  by  nearly  every  component  or  element 
through  changes  in  the  phase  and  amplitude  distribution  of 
the  beam.  Therefore,  if  realistic  estimates  of  the  per- 
formance of  candidate  laser  systems  are  to  be  obtained 
during  these  application  studies,  it  is  important  to  include 
the  effect  of  the  opt'cal  train  on  the  characteristics  of 
the  laser  beam. 

An  accurate  assessment  of  the  impact  of  the  ptical 
train  is  a tedious  and  difficult  calculation  requiring  a 
wave  optics  approach  and  sophisticated  analytical  tools  for 
modeling  each  of  the  elements  in  the  optical  train.  Obviously 
these  are  not  very  practical  for  use  in  systems  analysis 
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application  studies  in  which  a large  parameter  space  must 
be  investigated.  Therefore,  if  the  effect  of  the  optical 
train  is  to  be  included  in  such  analyses,  the  development 
of  simple  but  reasonably  accurate  models  is  required. 

One  such  model  has  been  developed  by  the  AFWL 
(Reference  5)  which  is  used  in  their  simplified  propagation 
code  (COMBO).  However,  its  use  is  somewhat  restrictive  in 
that  most  of  the  component  effects  are  left  for  the  user 
to  specify.  Prescriptions  for  scaling  these  effects  in 
terms  of  system  parameters  are  not  given.  In  addition,  the 
impact  of  the  optical  train  on  the  beam  quality  is  assessed 
in  terms  of  peak  intensity  reduction.  Based  on  work  under 
this  contract,  it  is  felt  that  an  assessment  in  terms  of 
the  far-field  power  distribution  provides  a better  charac- 
terization of  the  beam  quality  (see  Reference  1 and  Sec- 
tion 4.3  of  this  report). 

We  have  developed  a more  comprehensive  "simplified" 
optical  train  model  which  can  be  used  to  estimate  the 
degradation  in  system  performance  caused  by  the  optical 
components  required  to  direct  the  beam  from  the  laser 
device  to  the  transmitting  optic  (i.e.,  aerodynamic  window, 
mirrors,  beam  expander,  etc.).  The  intent  of  the  model  is 
not  to  provide  detailed  engineering  design  data  but  rather 
to  provide  the  systems  analyst  with  a rapid  assessment 
capability  which  will  allow  him  to  perform  more  realistic 
systems  analyses  of  candidate  HEL  systems.  Hence,  emphasis 
is  placed  on  formulating  simple  but  reasonably  accurate 
models  for  the  various  factors  influencing  the  performance 
of  the  laser  system.  The  approach  is  similar  to  that  taken 
by  the  AFWL  in  their  simplified  model  but  with  new  addi- 
tions and  modifications  to  make  the  model  more  flexible. 

The  model  has  been  implemented  in  a subroutine  called  OPTRAIN 
and  is  included  in  the  SAICOM  code  as  described  in  Section  3 
be! ow. 
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2.1  MODEL  FORMULATION 

Conceptually,  the  optical  train  may  be  viewed  as  a 
"black  box"  located  between  the  laser  device  and  the  trans- 
mitting optic.  The  input  to  the  box  is  the  output  beam 
from  the  laser.  The  model  operates  on  this  beam  according 
to  the  number  and  type  of  elements  in  the  train,  and  its 
output,  in  the  form  of  certain  performance  parameters, 
provides  the  input  to  an  atmospheric  propagation  model. 

The  performance  parameters  are:  (1)  the  power  available 

at  the  transmitting  optic,  (2)  beam  quality,  (3)  beam 
divergence,  and  (4)  beam  jitter.  In  addition,  the  actual 
diameter  of  the  beam  at  the  exit  aperture  is  computed  for 
use  in  the  propagation  model. 

The  essential  features  of  the  model  are  illustrated 
in  Figure  2.1.  The  optical  elements  and  components  include 
(1)  an  aerodynamic  window,  (2)  a beam  clipper,  (3)  several 
high  power  mirrors,  (4)  a beam  expander,  (5)  an  exit  aper- 
ture, and  (6)  a focusing  element  for  transmitting  the  beam 
to  the  target  plane.  In  most  practical  cases  the  beam 
expander  and  focusing  element  are  combined  into  a single 
telescope.  However,  it  is  convenient  for  the  present  anal- 
ysis to  separate  them  into  an  ideal  beam  expander  which 
simply  expands  the  beam  and  a separate  element  which  applies 
the  curvature  to  the  beam  phase  front  for  focusing  purposes. 

2.1.1  Input  Beam 

Although  it  is  not  an  optical  element  or  component, 
the  input  beam  is  necessarily  an  integral  part  of  the  model 
since  it  has  strong  influence  on  the  behavior  of  the  "est 
of  the  optical  train.  In  order  to  simplify  the  analysis 
of  cue  components  in  the  optical  train,  it  is  assumed  that 
the  beam  from  the  laser  is.  circular.  However,  because  of 
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the  widespread  use  of  unstable  resonator  configurations  to 
extract  power  from  high  power  devices,  the  possibility  of 
* central  obscuration  in  the  input  beam  1$  incorporated 
into  the  model,  the  inner  and  outer  diameters  of  the  input 
beam  are  cDg  and  Dg  respecti vely . An  unobscured  input  beam 
profile  can  be  realized  by  defining  e a 0.  Other  parameters 
required  to  specify  the  input  beam  are  (1)  the  power  in  the 
beam,  PL,  (2)  the  pulse  lengthi  At,  (3)  the  curvature  of 
the  phase  front,  R,  , and  the  beam  jitter,  0,. 

In  order  to  evaluate  the  effect  the  beam  has  on  the 
high  power  mirrors  in  the  optical  train,  it  is  also  nec- 
essary to  characterize  the  intensity  profile  entering  the 
optical  train.  This  is  achieved  by  specifying  both  the 
magnitude  and  scale  size  of  the  intensity  fluctuations, 
the  magnitude  being  specified  as  a fraction  of  the  average 
intensity  (AI/Iaye)  and  the  scale  size  as  a fraction  of 
the  beam  diameter  Uj/Dg).  Thus  the  input  beam  profile 
may  be  visualized  as  having  a uniform  intensity  Iaye  with 
fluctuations  AI  superimposed  upon  this  level. 

The  optical  quality  of  the  input  beam  is  specified 
by  a wavelength  scaling  factor,  nL>  Unfortunately,  this 
definition  of  beam  quality  does  not  allow  one  to  easily 
incorporate  the  other  beam  quality  degrading  factors  induced 
by  the  various  elements  of  the  optical  train.  Therefore 
this  parameter  is  converted  by  the  model  into  an  equivalent 
root-mean-square  phase  distortion,  a 

Two  alternative  approaches  for  this  conversion  were 
investigated  for  use  in  the  model.  The  first  approach  is 
to  define  the  equivalent  phase  distortion  so  as  to  match 
the  reduction  in  the  peak  intensity,  Ip.  According  to  the 
analysis  presented  in  Reference  6,  the  peak  intensity  reduc- 
tion is  (for  a gaussian  beam  containing  small  scale  random 
phase  variations) 
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where  Ig  is  the  ideal  on-axis  intensity.  According  to  the 
wavelength  scaling  definition  of  beam  quality 


so  that 

aL  * >/2tn(nL)  . v (1) 

It  should  be  noted  that  this  approach  is  equivalent  to  the 
power  scaling  definition  of  beam  quality  and  is  also  the 
approach  taken  by  the  AFWL  in  their  optical  system  model 
(Reference  5). 

The  second  approach  investigated  was  to  define  the 
equivalent  phase  distortion  to  match  the  power  in  the 
"bucket"  (RA/D  * 1)  predicted  by  wavelength  scaling.  This 
is  given  by  (Reference  1 ) : 

P ( RA/D)  * P0(l  - exp(-ir2/2nL2)}. 

Again,  following  the  analysis  of  Reference  4,  the  power  in 
the  bucket  for  a non-diffraction-limited  gaussion  beam  Is 
given  by 

P ( RA/D ) = PQ  exp(  -a2  ) (1  - exp  ( -tt2/2)  > 

assuming  ill  of  the  quality  degradation  can  be  characteri zed 
as  a scattering  loss.  Equating  these  expressions  and  solving 
for  yields, 


t 


1 - exp( -w */2n^) 
‘T-'txp'TTWy 


(?) 


Intuitively  it  would  seem  that  both  of  these  approaches 
should  yield  approximately  equivalent  results.  Quantita- 
tively however,  substantial  differences  were  found.  These 
are  illustrated  In  Figure  2.2  which  compares  the  phase  dis-  j 

tortion  computed  from  Equations  (1)  and  (2).  The  phase  J 

distortion  based  on  matching  the  peak  intensity  reduction  j 

is  always  considerably  higher  than  that  based  on  matching  j 

the  power-in-the-bucket,  e.g.,  over  a factor  of  2.5  for  a ; 

1.5  times  diffraction-limited  input  beam.  The  reasons  for 
these  differences  are  not  clear.  One  possible  explanation 
is  the  equivalence  of  wavelength  and  power  scaling  methods 

j 

of  defining  peak  intensity  reduction.  It  was  found  during 
previous  work  under  this  contract  (Reference  1)  that  power 
scaling  was  always  more  pessimistic  in  its  predictions  of 

the  beam  quality  than  wavelength  scaling.  Therefore,  ] 

matching  the  peak  Intensity  based  on  a method  that  is 

equivalent  to  power  scaling  might  be  expected  to  produce 

more  pessimistic  results  (i.e.,  higher  crL)  than  matching 

the  oower-i n-the-bucket  based  on  wavelength  scaling. 


Whatever  the  reasons  for  the  differences,  we 
recommend  the  use  of  Equation  (2)  because  the  power-in-the- 
bucket  is  a more  meaningful  measure  of  laser  performance 
than  peak  intensity  reduction.  The  latter  suffers  from 
its  sensitivity  to  defocusing  errors  and  difficulty  in 
experimental  measurement. 

2.1.2  Aerodynamic  Window 

The  first  component  in  the  optical  train  is  the 
aerojynamic  window  isolating  the  optical  cavity  of  the  laser 
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from  the  atmosphere.  Basically  there  ere  two  types  of 
aerodynamic  windows  currently  being  considered  for  Hit 
applications:  (1)  focused  end  (2)  unfocused.  In  both 

cases,  If  the  window  Is  properly  designed,  there  should 


be  little  or  no  power  loss  when  the  beam  passes  through 
It.  In  eddltlon,  it  Is  also  assumed  that  no  spherical 
distortion  to  the  phase  front  will  be  Induced  by  the  aero- 
dynamic .-Indow.  The  shock  waves  supporting  the  pressure 
rise  and  the  turbulence  generated  by  the  aerodynamic  win- 
dow will,  however,  Induce  a loss  In  the  beam  quality.  An 
expression  for  this  loss  Is  given  by  {Reference  4) 

0 2 

A 1/ 1 ■ 2.P4  x 10-10 

where  Ap  is  density  change  across  the  shear  laye*  and  p& 

is  a reference  density  at  STP.  Again  this  form*;! Ism  of  1 

beam  quality  degradation  Is  not  convenient  for  use  In  the 
model  and  is  converted  into  an  equivalent  phase  distortion, 
o£w  by  use  of  the  Strehl  formula  (Reference  4): 

M/I  = j 

i 

i 

j 

Solving  for  the  phase  distortion  yields  j 

'«  ■ >-69  * 10-5  (t  f;)-  ] 

1 

j 

A typical  value  of  for  a 10  cm  beam  from  a chemical 

a n 

laser  (X  - 3.8  ym),  assuming  up/ps  = 0.25  for  the  density 

variations,  is  baw  * 0.11  radians.  This  corresponds  to  an 

intensity  reduction  of  AI/I  5:  i.2S,  I 


For  the  focused  aerodynamic  window,  the  beam 
diameter  will  typically  be  very  small  {<0.1  cm)  so  that 
the  phase  distortion  can  be  neglected,  l.e.,  <?4W  « 0. 

2.1.3  Seam  Clipper 

For  MEL  systems  in  which  the  laser  device  and  the 
pointer/tracker  are  not  closely  coupled,  a beam  clipper  or 
scraper  will  probably  be  employed  near  the  pointer/tracker 
to  control  the  site  of  the  high  power  beam  entering  the 
transmitter.  The  beam  clipper  will  influence  both  the 
power  available  at  the  transmitting  optic  and  the  bean 
quality,  the  latter  by  limiting  the  maximum  diameter  at 
the  transmitting  optic. 

The  loss  in  power  at  the  beam  clipper  is  due  pri- 
marily to  diffractive  spreading  of  the  beam  as  It  propagates 
between  the  aerodynamic  window  and  the  beam  clipper.  In 
genera'.,  the  diffractive  spreading  will  depend  upon  the 
phase  and  amplitude  aberrations  in  the  input  beam  profile, 

An  accurate  assessment  of  the  actual  losses  requires  a 
detailed  calculation  involving  numerical  solutions  of  the 
Fresnel  diffraction  integral  or  spatial  frequency  approaches 
to  ootical  propagation  (l.e.,  Fourier  Optics).  These  tech- 
niques are  very  time  consuming  and  wero  net  considered  as 
viable  candidates  for  the  "simplified"  model.  Instead,  an 
empirical  approach  based  on  curve  fits  of  available  data 
in  the  literature  was  taken.  Unfortunately  very  little 
data  wore  available  for  examination.  Only  two  sources  were 
found  (References  4 and  7).  The  data  obtained  from  these 
sources  are  illustrated  in  Figure  2.3  which  shows  the  vari- 
ation in  power  diffracted  into  the  geometric  shadow  of  a 
uniformly  illuminated  aperture  as  a function  of  the  propaga- 
tion distance  from  the  aperture.  The  AFWl  data  were  calcula 
ted  for  a one-dimensional  beam  profile  while  the  MAC  (Hughes 


Figure  2.3.  Diffractive  Spillover  Losses 


Aircraft  Company)  calculations  were  performed  for  a found 
beam.  Therefore , the  data  must  be  adjusted  for  these 
differences  before  a direct  comparison  can  be  made. 

Intuitively  one  would  expect  a factor  of  two 
difference  between  the  data  based  on  simple  geometrical 
arguments,  for  example,  a square  beam  profile  would  have 
twice  the  loss  of  a one-dimensional  profile  because  of  the 
two  additional  edges,  for  a round  beam,  however,  the  ratio 
is  not  quite  as  straightforward.  Assuming  circular  symmetry, 
the  fractional  power  loss  is 


where  a is  the  radius  of  the  aperture  and  l0  is  the  initial 
Intensity.  The  major  contribution  to  the  integral  occurs 
near  the  edge  of  the  aperture.  Thus  over  this  region  the 
radius  can  be  taken  as  being  approximately  constant  {•  a) 
and  removed  from  the  integral  without  affecting  its  value, 
i ,e. , 


‘2-0  * aT 


«• 

- f Ho 

Q J 


) da . 


For  the  one-dimensional  beam  profile  the  fractional  power 
loss  is  given  by 


J Kx 


) dx  (1) 


kl-D 
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to  that  If  one  further  assumes  that  the  Intensity  prof  Hot 
• t tit*  aperture  plane  art  the  tame  in  either  case,  i.e., 
!{r)  * l(x),  then  one  again  gets  3 factor  of  two  difference 
between  the  potier  losses. 

The  actual  ratio  indicated  by  the  data  was  somewhat 
higher  than  a factor  of  <t  as  shown  by  the  corresponding 
curve  fits  of  the  data,  toth  sets  cf  data  correlated  very 
well  as  a function  of  the  reciprocal  of  the  square  root  of 
the  Fresnel  number.  However,  the  proportionality  constant 
for  the  HAC  data  was  found  to  be  23.6  versus  li.l  for  the 
AFWl  data. 

In  order  to  use  these  results  in  the  optica)  train 
mode),  the  power  loss  was  assumed  to  be  related  to  an 
effective  increase  in  the  beam  diameter  at  the  aperture 
plane,  i.e., 

n - i.,.)1'*  i - 

*0  '"T 

where  0B  is  the  new  diameter  of  the  beam  at  the  beam 
B°  2 

clipper  station  and  FN  « Dj/4Alc  is  the  propagation  Fresnel 
number.  It  was  also  assumed,  for  the  case  of  an  obscured 
beam,  that  the  hole  in  the  beam  filled  in  as  much  as  the 
diameter  of  the  hale  would  grow  were  it  a uniformly  illu- 
minated aperture  (i.e.,  habinet's  Principle).  Thus  the 
inner  diameter  of  the  obscured  beam,  0^^ , is  computed  by 

°Bi  * t0B  <l  ' 0.112/cV^j). 


If  &c  Is  tilt  diameter  of  the  beam  clipper,  then  the  power 
transmission  through  the  aperture  Is 


T 


( 1*0  If  o,#  s oe 

j <VV  U \ * °‘ 
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In  the  litter  cue*  the  beam  diameter  Is  set  equal  to  the 
clipper  diameter  for  the  remaining  calculations  In  the 
model . 

2.1.4  Hlqh  Power  Mirrors 

The  mirrors  In  the  optical  train  will  affect  all 
the  performance  parameters.  Although  It  Is  not  necessary 
we  assume  that  all  of  the  mirrors  are  Identical  In  order 
to  reduce  the  number  of  Inputs. 

A mirror  alters  the  amplitude  of  the  beam  by 
absorbing  a fraction  of  the  Incident  radiation.  For  a 
series  of  N mirrors  each  having  a surface  reflectivity  A, 
the  transmission  factor  Is 

tm  . 

It  should  be  noted  that  In  the  present  analysis  N Includes 
alt  of  the  mirrors  In  the  optical  train,  l.e.,  relay, 
steering,  secondary  and  primary. 

The  beam  quality  degradation  Is  due  to  (1)  surface 
roughness  and  manufacturing  errors  In  the  mirror  figure, 
and  (2)  surface  distortions  due  to  absorption  of  power  from 
the  Incident  beam.  The  former  are  Independent  of  the  power 
In  the  buam  while  the  latter  Is  a function  of  both  the 
power  and  the  Irradlance  distribution.  In  addition.  If 
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cooled  mirrors  are  used,  there  will  also  be  distortions 
induced  by  the  pressure  variations  inside  the  coolant 
passages.  However,  these  distortions  will  depend  upon  the 
details  of  the  mirror  construction  and  their  characteriza- 
tion is  felt  to  be  beyond  the  scope  of  this  model.  There- 
fore, for  this  analysis,  distortions  of  this  type  will  be 
included  in  (1)  above. 


The  beam  quality  degradation  induced  by  the  mirrors 
is  incorporated  into  the  optical  train  model  through  an 
equivalent  phase  variance,  c?M  where 


a2  = Naf  + - D al  + Cj/M2]2. 

The  first  term  on  the  right  hand  side  of  the  above  expres- 
sion represents  the  uncorrel ated  sum  of  the  phase  distortion 
due  to  fabrication  errors  such  as  surface  roughness,  figure 
error,  coolant  passage  distortion,  etc.  The  remaining  terms 
represent  the  correl ated  sum  of  the  irradiance  mapping  phase 
distortion.  The  distortion  of  the  primary  mirror  is  reduced 
by  the  square  of  the  magnification  of  the  beam  expander  to 
account  for  its  reduced  thermal  loading. 


In  practice,  is  not  known  but  is  instead  specified 
as  a tolerance  on  the  manufacturing  process.  Therefore, 
for  modeling  purposes  this  parameter  was  left  to  the  user 
to  specify. 

The  approach  taken  to  model  Oj  was  to  assume  that 
the  distortions  scale  directly  with  the  irradiance  fluctua- 
tions. That  is, 

oj  ='  KI 


rms 


where  K is  a constant  which  depends  upon  the  physical 
parameters  of  the  mirror,  the  cooling  scheme,  etc.  If  we 
assume  that  the  peak  to  peak  Intensity  fluctuations  speci- 
fied in  the  Input  beam  description  are  uniformly  distributed 
and  random  In  character,  then 


I 

rms 
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where,  for  a circular  beam, 
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Expressions  for  the  evaluation  of  K depend  upon 
whether  the  mirror  is  water  cooled  or  not,  i.e.,  for  cooled 
optics 

K(m2/watt)  = 5 x 10"14  (1  - R)  (2tt/X) 


whereas,  for  uncooled  optics 


where 


K ■ 4ir  (1  - R)  aAt/XCp 


a = the  thermal  expansion  coefficient  of  the 
mirror  material 

C = the  specific  heat  of  the  mirror  material 


The  expression  for  cooled  optics  is  based  on  the 
NPT/Chemical  Laser  Compatibility  Study  conducted  by 
Hughes  Aircraft  Company  for  NSWC  (Reference  8)  in  which 
similar  expressions  were  employed  to  compute  the  phase 
distortion  caused  by  the  NACL  irradlance  profile.  The 
expression  for  uncooled  optics  is  based  on  a simplified 
one-dimensional  heat  transfer  analysis  (Reference  9). 

Local  surface  distortion  is  not  the  only  performance 
degrading  factor  caused  by  the  finite  absorptivity  of  the 
mirrors.  A bending  distortion  is  also  induced  by  the 
differential  growth  of  the  front  and  rear  mirror  surfaces. 
For  properly  designed,  cooled  mirror  configurations,  this 
distortion  can  be  kept  small.  However,  for  uncooled  mirrors 
it  may  be  important. 

To  first  order,  the  bending  distortion  is  primarily 
a function  of  the  total  beam  power  and  not  the  irradiance 
distribution.  In  addition  it  produces  mostly  a spherical 
phase  front  distortion.  Thus  in  the  optical  train  model, 
we  compute  this  distortion  mode  as  a beam  divergence  instead 
of  a quality  loss  since,  in  theory,  it  could  be  corrected  by 
the  focusing  optic  (if  detected). 

For  an  uncooled  mirror  that  was  initially  flat,  the 
change  in  the  focal  length  with  time  can  be  approximated 
as  (see  Reference  9) 

_ p|  * COS  a /DgK  \1/2 

T " 2.44  PL  (1  - R)  a \ At  / 


where 


i is  the  mirror  thickness 
0 is  the  beam  ang'ie  of  incidence  on  mirror 


K Is  the  thermal  conductivity  of  the  mirror 
material 

f is  the  focal  length  of  the  distorted  mirror. 

The  mechanism  by  which  the  defocusing  errors 
induced  by  the  individual  mirrors  are  accumulated  through- 
out the  optical  train  is  explained  later. 

2.1.5  Beam  Expander 

In  the  formulation  of  the  beam  expander  model,  we 
allow  for  two  types:  (1)  an  on-axis  or  (2)  an  off-axis 

system.  Either  type  system  affects  the  available  power 
at  the  transmitter,  the  beam  quality,  beam  jitter  and  beam 
divergence.  A discussion  of  the  methodology  for  computing 
the  transmission  factor  and  beam  quality  is  given  below. 
The  beam  jitter  and  divergence  are  discussed  in  a later 
section. 

The  power  transmission  through  the  beam  expander 
is  computed  by  projecting  the  exit  aperture  (defined  by 
the  user)  onto  the  plane  of  the  beam  clipper  station. 
Simple  geometry  then  allows  one  to  compute  the  power  lost 
due  to  a mismatch  between  the  beam  diameters  and  exit 
aperture  diameters.  The  formulas  are  illustrated  in 
Table  2.1  for  both  the  on-axis  and  off-axis  cases. 

The  beam  quality  is  not  affected  by  the  off-axis 
system  except  to  the  extent  that  the  limiting  diameter  is 
MDg  < for  spot  size  calculations  at  the  target  plane. 

This  is  also  true  of  the  on-axis  case.  However,  the  on- 
axis  system  also  introduces  a central  obscuration  e,  which 
reduces  the  beam  quality.  The  expressions  used  in  the 
model  to  compute  the  final  obscuration  of  the  beam  leaving 
the  optical  train  are  also  given  in  Table  2.1. 
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The  effect  of  the  obscuration  on  the  beam  quality 
is  again  computed  by  relating  i*  to  an  equivalent  phase 
distortion.  This  relationship  was  derived  empirically 
based  on  the  results  of  the  previous  work  under  this  con- 
tract (Reference  2),  and  is  illustrated  graphically  in 
Figure  2.4.  Within  the  computer  code,  this  relationship  Is 
represented  by  a third  order  polynomial  developed  from  a 
"least-squares"  regression  analysis  of  the  data  presented 
in  Figure  2.4.  Also  shown  on  the  figure  is  the  effective 
scale  size  of  the  phase  distortion.  Its  use  in  the  model 
will  be  explained  below. 


2,1.6  Exit  Aperture 

Provisions  were  also  made  in  the  optical  train 
model  for  an  exit  aperture  downbeam  of  the  pointer/tracker. 
It  can  be  either  a material  window  or  open  port,  with  an 
aerodynamic  curtain  protecting  the  optics  from  the  environ- 
ment, based  on  user  specified  option  parameters. 

In  either  case  the  loss  in  power  through  this  exit 
aperture  will  be 


where 

$ = absorption  coefficient  of  material  window 

lw  « thickness  of  material  window 

A « area  of  struts  supporting  window  or  secondary 
mirror  in  the  case  of  an  on-axis  beam  expander 

2 2 

Ag  = beam  area  (ttM  Dg  /4). 
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Note  that  an  open  port  can  be  simulated  by  setting  8 ■ 0. 

For  simplicity*  the  phase  distortion  Induced  by 
the  material  window  Is  assumed  to  be  proportional  to  the 
local  Intensity  fluctuations  In  the  beam  profile.  That  Is 


■ (t) 


S&tl, 


i)  « ♦ If 
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where  (in  addition  to  those  parameters  previously  defined) 
n is  the  refractive  Index  and  T is  the  temperature.  The 
f jnlfication  appears  because  of  the  assumed  location  of 
the  window  downstream  of  the  beam  expander.  In  order  to 
minimize  the  number  of  inputs  to  the  model,  the  raterlal 
constants  were  lumped  into  a single  parameter  defined  as 


(n  - 1)  a + 


12  3 

whu  varies  between  4 x 10“  m /j  for  the  fluoride  win- 
dows af2,  MgF2,  SrF2  ) to  15  x 10"12  m3/j  for  the  salt 
windows  (NaCl , KC1).  A value  of  5 x 10"12  m3/j  was  hard- 
wired ii  .o  the  model  as  being  representative  of  current 
window  technology. 

The  presence  of  struts  in  the  beampath  will  also 
cause  a reduction  in  the  beam  quality.  To  first  order, 
this  loss  will  be  directly  proportional  to  the  area  of  the 
beam  blocked  by  the  struts.  Consider,  for  example,  the 
on-axis  intensity  which,  in  the  context  of  scalar  diffrac- 
tion theory,  is  given  by  th-:  Integral  of  the  complex  field 
leaving  the  transmitting  aperture,  i.e., 


I(o) 


l ff 


(x,y)  d x dy 
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Thus,  for  a constant  field  aperture  strength,  the  integral 
Is  simply  the  area  of  the  clear  aperture.  That  is 


Ko)  ' (Aj  - A,)* 


so  that  the  intensity  reduction  relative  to  no  struts  is 
simply 


For  modeling  purposes,  this  intensity  reduction  was  con- 
sidered to  be  separable  into  (1)  a power  loss  due  to  the 
transmission  of  the  aperture,  Ts,  and  (2)  a beam  quality 
loss  characterized  by  wide  angle  scattering,  (1  - A^/Ag). 
The  latter  can  be  related  to  an  equivalent  phase  distortion 
via  the  Strehl  equation 


°s  * (As/Ag)**  . 


The  phase  distortion  induced  by  an  open  port  is 
not  intensity  dependent  but  instead  a function  of  the  tur 
bulence  level  inside  the  beam  expander  and  aerodynamic 
curtain.  Because  of  the  presence  of  turbulence,  this 
effect  is  treated  as  an  adoitional  source  of  beam  jitter 
in  the  model.  The  magnitude  of  the  jitter  is  estimated 
from  (Reference  8) 


Oj  * 2.14  x 10’ 


Le  AT6 
X D* 
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where  Lf  Is  the  path  length  through  the  bean  expander  and 
AT  Is  the  magnitude  of  the  temperature  fluctuations  within 
the  beam  expander.  Some  representative  values  of  the 
jitter  are  Illustrated  In  Figure  2.5  for  a 0,?  meter  trans- 
mitter diameter.  For  example,  a 1°C  temperature  fluctuation 
produces  approximately  4.S  yrad  of  beam  jitter  for  a 
3.8  pm  wavelength  beam. 

2.2  ASSESSMENT  OF  OVERALL  PERFORMANCE 

The  previous  discussion  was  primarily  concerned 
with  the  Influence  each  of  the  elements  has  on  the 
characteristics  of  the  laser  beam  as  It  traverses  the 
optical  train.  In  the  following  section  we  show  how  these 
Individual  effects  are  accumulated  to  arrive  at  an  overall 
assessment  of  system  performance. 

2.2,1  Transmitted  Power 

The.  power  available  at  the  transmitting  optic  Is 
simply  the  input  power  from  the  laser  device  modified  to 
account  for  all  of  the  power  losses  that  have  occurred 
along  the  optical  train.  That  is 


■ T \ 


where  T Is  an  overall  transmission  factor.  Since  the 
power  losses  are  multiplicative,  the  overall  transmission 
factor  Is 


T * <V  <V  (Tbe)  (Te) 
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where  T£  is  the  transmission  factor  of  the  beam  clipper. 

It  the  transmission  factor  of  til  th#  high  power  mir- 
rors, TQ£  Is  the  transmission  ftctor  of  the  beam  expander* 
and  Is  the  transmission  ftctor  of  the  exit  aperture. 
Formulas  for  til  of  these  factors  were  given  previously. 

2.2.2  Beam  Quality 

The  quality  of  the  bean  at  the  transmitting  optic 
Is  characterised  by  the  two  beam  quality  parameters*  mj 
and  m2  (see  Reference  l and  Section  4.3  of  this  report). 
Unfortunately*  the  calculation  of  m^  and  m2  Is  not  simple. 

The  basic  problem  Is  to  combine  the  effect  of  all  the 
quality  degrading  factors  occurring  within  the  optical 
train.  This  Is  somewhat  complicated  because  the  phase 
distortions  induced  by  the  individual  components  do  not 
add  up  in  a straightforward  manner  but  instead  depend*  in 
a complicated  manner,  upon  the  distribution  as  well  as  the 
magnitude  of  the  phase  distortion.  For  example,  a "smooth" 

\ distortion  behaves  differently  than  a very  "rough"  distor- 

tion even  though  their  magnitudes  are  the  same. 

To  circumvent  this  problem,  a statistical  approach 
Is  used  within  the  model  to  accumulate  the  phase  distortions. 
This  approach  is  somewhat  loosely  based  on  the  theoretical 
Investigations  of  nondi ffractlon-1 1ml ted  gaussian  beams  by 
B,  Hogge  at  the  AFWL.  Briefly,  the  approach  is  to  consider 
the  overall  beam  profile  at  the  focal  plane  as  being  com- 
posed of  two  gaussian  beam  profiles  of  different  relative 
amplitudes  and  widths.  That  Is,  a certain  fraction  of  the 
energy  in  the  beam  is  propagated  completely  unaffected  by 
the  phase  distortion,  l.e.. 


i 

i 
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where  Wf  It  the  diffraction  11t»Ued  waist  parameter,  tQ 
Is  tht  on-axis  Intensity  and  o*  It  tht  variance  of  tht 
phtst  distortion  (assumed  to  « geussian  random  variable). 
Tht  remaining  tntrgy  Is  smeared  or  spread  by  tht  phtst 
distortion  into  t somewhat  larger  bttti  profile,  i.e.. 


where  © is  tht  angular  sprttd  of  tht  scttttrtd  beam. 

Tht  problem  thus  reducts  to  chtrtcterit ing  the 
overall  phtse  variance,  a,  tnd  beam  spread  parameters,  ©, 
that  "best*  model  the  summation  of  all  of  the  phase 
distortions  in  the  optical  train, 

A true  characterization  of  o and  © is,  in  reality, 
not  possible.  Some  theoretical  arguments  can  te  made  for 
Isolated  phase  ditributions  such  as  the  random  gaussian 
noise  model  Investigated  by  Hogge,  In  the  real  world,  how- 
ever, optical  systems  do  not  produce  such  amiable  distor- 
tions. Nevertheless,  using  these  results  as  being  at  least 
qualitatively  correct,  we  compute  the  total  phase  distortion 
as  simply  the  uncorrected  sum  of  each  of  the  individual 
distortions,  i.e.. 


Instead  of  computing  the  beam  spread,  o, 
convenient  tc  compute  an  effective  scale  size  or 
length  of  the  overall  phase  distortion,  t.  This 
is  related  to  the  beam  spread  via  0*  « 


It  Is  more 
correlation 
scale  size 
In  the 
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aforementioned  work  of  Nogge,  It  was  noted  that  the 
effective  correlation  length  for  a number  of  sources  of 
phase  distortion  was  found  to  be  an  average  of  the  correla- 
tion lengths  of  the  Individual  sources,  each  weighted  by 
Its  respective  variance.  Again,  taking  the  Inductive  leap, 
we  compute  L from 

1 * [k®?  * *AW°AW  4 lNf1H  4 tB£°8E  ' *$®S  4 *W°w]  ’ 

v 

In  prartlce  very  little  Is  known  about  the  magnitude  of 
the  scale  sixes  characterising  the  Individual  components. 
Therefore,  for  the  present  model,  *j  have  made  some  addi- 
tional assumptions  regarding  their  respective  sixes.  For 
example,  the  Irradlance  mapping  phase  distortions  such  as 
o | and  are  assumed  to  have  the  same  scale  size  as  the 
intensity  fluctuations  in  the  Input  beam  profile.  For 
other  components,  such  as  the  laser  device,  the  aerodyanmlc 
window  and  the  struts,  we  assume  that  the  scale  size  is 
zero.  In  effect  this  is  assuming  that  each  of  these  com- 
ponents scatters  the  energy  In  the  beam  beyond  a usable 
radius  In  the  focal  plane  and  therefore  is  somewhat  conserva- 
tive. The  scale  size  of  the  beam  expander  was  determined 
empirically  since  the  effect  of  obscuration  on  the  near- 
field  beam  profile  was  readily  calculated  (see  Figure  2.4), 
Within  the  model  this  relationship  is  represented  by  a 
third  order  polynomial. 

The  computed  values  of  the  phase  variance  and 
scale  size  are  used  to  evaluate  the  power  distribution  at 
two  radial  positions  in  the  far-fleld,  namely  r « w^  and 
fg  » 2wf.  These  power  points  are  then  used  to  generate 
the  beam  lity  parameters  aij  and  m.,.  The  procedure  for 
doing  this  Is  outlined  in  Reference  1. 

2**8 
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2.2.3  Beam  Jitter 

The  beam  jitter  is  accumulated  throughout  the 

optical  train  by  assuming  that  the  sources  of  beam  jitter 

are  independent  of  each  other.  Thus  they  may  be  root-sum- 

squared  to  gee  the  total  jitter.  The  primary  consideration 

for  the  beam  jitter  calculation  is  whether  or  not  the 

source  of  the  jitter  is  upbeam  or  downbeam  of  the  beam 

expander.  The  beam  expander  reduces  the  jitter  by  a 
2 

factor  of  1/M  where  M is  the  magnification. 

The  sources  of  beam  jitter  considered  by  the  model 


are 


Tracking  jitter, OrR 

Boresight  jitter,  GgS 

Autoalignment  system  jitter,  oAA 

Servo  jitter,  0S 

Device  jitter,  0L 

Exit  aperture  induced  jitter,  . 

The  jitter  of  the  input  beam  and  the  jitter  induced  by  the 
beam  steering  mirror  in  the  autoalignment  system  are 
assumed  to  occur  before  the  beam  expander.  The  remaining 
jitter  sources  are  assumed  to  occur  after  the  beam  expander. 
Thus  the  total  jitter  leaving  the  optical  train  is  given  by: 


>0  * <®TR  + ®BS  + ®S  + ®E  + <®L  + ®AA>/m2>,S- 


2.2.4  Beam  Divergence 

The  beam  divergence  represents  the  spherical  phase 
curvature  in  the  beam  as  it  enters  the  transmitting  optic 
and  results  in  a larger  spot  size  at  the  target  than  one 
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would  expect  had  the  beam  been  collimated  when  It  entered 
the  transmitting  optic.  As  an  example,  consider  a colli- 
mated ideal  gausslan  beam  which  has  been  focused  to  pro- 
duce  a spot  size  of  wf  (actually  a radius  to  the  e"c 
intensity  point)  at  a target.  If  the  beam  is  instead 
diverging  (or  converging)  with  a radius  of  curvature  -Rg 
(+Rg)  as  it  is  reflected  from  the  same  transmitting  optic, 
the  new  spot  size  w|  is  larger  than  the  original  by 


2 

where  Dg  is  the  (e  ) diameter  of  the  ideal  gaussian  beam 
and  X is  the  wavelength.  For  example,  0.1X  spherical  phase 
error  will  increase  the  spot  size  approximately  20%  for  a 
100  cm  diameter  optic  transmitting  at  X = 3.8  y.  For  the 
present  study  we  assume  that  any  phase  curvature  in  the 
output  beam  from  the  optical  train  behaves  in  a similar 
manner. 

The  change  in  the  curvature  of  the  phase  front  of 
the  beam  as  it  propagates  through  the  optical  train  Is 
computed  by  employing  a simplified  geometric  optics  calcula- 
tion through  the  optical  elements  preceding  the  beam 
expander.  The  usual  simplifying  assumptions,  i.e.,  paraxial 
rays,  thick  lenses,  etc.,  are  made  in  order  to  make  the 
calculation  tractable. 

Following  the  ray  matrix  approach  of  Reference  10 
we  first  obtain  the  equivalent  ray  matrix  of  the  optical 
train  illustrated  in  Figure  2.6  by  multiplying  together  the 
individual  ray  matrices  of  the  individual  elements,  I.e., 
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where  (going  from  right  to  left)  the  first  matrix  represents 
the  propagation  distance  between  the  aerodynamic  window  and 
the  beam  clipper,  the  second  matrix  represents  the  reflec- 
tion of  the  beam  from  the  first  relay  mirror  In  the  optical 
train,  the  third  matrix  represents  the  propagation  distance 
between  the  first  relay  mirror  and  the  second  relay  mirror, 
the  fourth  matrix  represents  the  reflection  of  the  beam 
from  the  second  relay  mirror,  and  so  on  until  the  beam  Is 
incident  upon  the  secondary  mirror.  At  this  point  we 
assume  that  the  beam  expander  simply  e>pands  the  beam  by 
the  specified  magnification.  The  cur  .ture  of  the  phase 
front  at  this  point,  Rg,  is  related  to  the  phase  front 
curvature  of  the  input  beam,  RL,  by 


Rb  = M 


ARl  + B 
CRl  + D 


The  computed  value  of  Rg  is  then  used  to  compute  the  new 
spot  size  of  the  focused  beam.  In  order  to  transfer  this 
information  to  the  propagation  model,  the  beam  quality 
parameter,  m2,  which  characterizes  the  spread  of  the  beam 
profile  in  the  abs°nce  of  other  effects,  such  as  thermal 
blooming,  turbulence  and  beam  jitter,  is  internally  adjusted 
by  the  model  to  reflect  the  increased  beam  spread,  i.e., 


2.3  SAMPLE  CALCULATIONS 

As  an  Illustration  of  how  the  optical  train 
influences  the  performance  of  a laser  system,  we  present 
below  a sample  calculation  for  the  optical  train  schemat- 
ically Illustrated  In  Figure  2.7. 

The  laser  device  Is  assumed  to  be  operating  In  an 
unstable  confocal  resonator  mode  with  the  power  coupled 
out  of  the  device  via  a scraper  mirror  located  in  front  of 
the  convex  optic  of  the  resonator.  Accordingly,  the  output 
beam  contains  a central  obscuration.  The  parameters 
assumed  for  this  beam  are: 

Beam  Diamter  ■ 0.1  meters 

Power  = 400  kw 

Wavelength  = 3.8  microns 

Obscuration  = 0.5 

In  addition,  it  is  assumed  that  the  beam  is  essentially 
collimated  ( = »)  and  contains  peak-to-peak  intensity 
fluctuations  that  are  50  percent  of  the  average  intensity 
(AI/Iave  = 0.5),  The  scale  size  of  these  fluctuations  is 
taken  to  be  1 cm  or  tj/Dg  = 0.1.  We  also  assume  that  the 
quality  of  the  beam  leaving  the  device  is  characterized  as 
being  1.2  times  diffraction-limited  (nL  = 1.2).  This 
corresponds  to  an  equivalent  rms  phase  distortion  in  the 
beam  of  oL  = 0.16  radians  or  about  A/40  (see  Figure  2.2). 

The  output  beam  from  the  laser  device  is  focused 
through  an  aerodynamic  window  to  a collimating  mirror. 

Since  the  beam  size  is  very  small  (typically  < 0.1  cm)  at 
the  aerodynamic  window,  we  neglect  any  phase  distortion 
Induced  by  the  aerodynamic  window  (i.e.,  oaw  * 0). 
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For  the  exempte»  we  assume  that  the  beam  expander 
Is  located  SO  meters  from  the  aerodynamic  window.  Over 
this  distance,  diffraction  effects  cause  the  beam  to 
expand  somewhat.  The  Fresnel  number  for  this  propagation 
distance  is 


«! 


13.16 


so  that  the  beam  diameters  at  the  entrance  to  the  beam 
expander  are 


!B0  " 1“-  0.112//F'N 


= 0 . 1 03  m 


and 


L 0.112\ 

r'^r) 


Og  * 0.047  m, 


In  this  example,  we  have  assumed  that  the  diameter  of  the 
beam  clipper  is  16  cm  so  that  no  power  loss  due  to  clipping 
occurs , 1 .e. , Tc  * 1.0. 

The  nine  (9)  mirrors  employed  in  the  optical  train 
are  assumed  to  be  water  cooled  with  a reflectivity  of 
* 0.986.  We  also  take  the  fabrication  error  as  being 
Xy/8  where  Xy  is  0.564  microns.  Accordingly,  the  power 
lost  to  the  beam  due  to  the  finite  absorptivity  of  the 
mirrors  is 

Tm  = (0.986)9  = 0.88. 
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The  phase  distortion  Induced  by  each  mirror  depends  upon 
the  rms  Intensity  fluctuations,  which  for  the  case  being 
considered  are 

« 7.35  x 106  w/m2 

12 

so  that 

a.  « 5 x 10"14  {1  - R ) ( 2tt/X ) I * 0.0085  radians.  j 

* 1 

3 

i 

-j 

Combining  this  with  the  fabrication  error  yields  (for  all 
of  the  mirrors  In  the  optical  train) 


4PL 

rms  ^ 


°M  * 9of  * 


[8<?j  + Oj/M2]2  « 0.3454  radians. 


The  beam  expander  is  an  on-axis  system  of  magnifica- 
tion. M:  - 4 . 3 " i . The  inner  and  outer  diameters  of  the 

clear  output  aperture  are  0+  3 0.7  and  0.  3 0.252  meters 

zo  t1 

respectively.  Since  the  beam  entering  the  expander  is  not 

perfectly  matched  to  the  reduced  clear  aperture  dimensions 

(i.e.,  D*  = 0.16  and  D*  3 0.576),  there  is  a power  less 
*o  *1 

induced  by  the  beam  ****:ander.  From  the  geometry  of  the 
case  being  considered,  ohis  loss  is  computed  to  be  approxi- 
mately 13  percent  or  TB£  ^ 0,87.  In  addition,  tne  beam 
expander  has  altered  the  obscuration  of  the  input  beam 
because  of  clipping.  final  obscuration  of  the  beam 

leaving  the  beam  expa  i :ar  is  e = 0.0576/10.3  3 0.558.  Note 
that  the  final  obscuration  is  not  the  observation  of  the 


telescope  (Dt  /Dt  « 0,36)  because  the  Input  bean  Mas  not 
*1 

large  enough  to  fill  the  transmitting  optic.  The  effective 
phase  distortion  of  the  actual  bean  Is  therefore  much 
larger  than  one  would  compute  based  on  the  obscuration  of 
the  telescope.  For  the  example  being  considered,  this 
equivalent  phase  distortion  Is  o&E  « 0.72  radians. 

The  exit  aperture  was  assumed  to  be  an  open  port 
so  that  the  only  power  loss  is  that  due  to  blockage  by 
the  struts  required  to  support  the  secondary  mirror  of  the 
beam  expander.  Thus,  T£  ■ 0.9S.  The  quality  loss  due  to 
the  struts  Is 

\ 

°s  " ^ * 0.2236  radians.  :~ 

Within  the  model,  the  open  port  would  also  Induce  jitter 
Into  the  beam.  However,  for  the  purposes  of  this  example, 
we  have  Ignored  any  performance  loss  caused  by  beam  jitter. 

Combining  the  above  losses,  we  get  for  the  overall 
transmission  factor 

T - Tc  * TBE  * TE  - °-727- 

Thus,  the  power  from  the  device  that  is  available  for 
propagation  to  "the  target  is 

/ 

Pt  - 0.727  PL  « 290  kw. 

The  tots!  distortion  Is 

a * (°m  + °BE  + °S  + 0«*33\  *1.06  rads. 


j 
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This  phase  distortion  Is  characterized  by  a scale  size 

w*  ^ °«  + 4 l(N ' 1,01  * ’ 0ll9Sl  * 0-114 

In  torsts  of  tht  two  parameter  descriptions  of  the  been  qual- 
ity employed  by  the  model  m^  » 0.3536  end  «2  • 1,.  1076.  For 
this  example,  there  were  no  defocusing  errors  Induced  by 
any  of  the  components  In  the  optical  train  so  that  no  adjust- 
ment was  made  to  the  beam  spread  parameter  m^.  Note  that  the 
constraints  In  the  last  two  equations  (0.433  and  0.196)  are 

needed  to  ensure  the  proper  conversion  of  a and  k-  to  m.  and 

m B 

In  order  to  quantify  the  Impact  of  the  optical 

tra1nt  we  evaluate  the  power  delivered  to  a particular  spot 

In  the  focal  plane,  namely  r*  ■»  fX/D.  . In  the  absence  of 

*o 

atmospheric  effects  and  beam  jitter,  this  power  Is  simply 


P(r*) 


tml  | 1 


*2  /°b\ 
* !T  I 5 7*1 
m2  \ J 


where  Dg  Is  the  actual  diameter  of  the  beam  on  the  trans- 
mitting optic,  l.e.,  0B  » N Dft  < D«.  . 

0 B0  lo 

The  results  of  this  calculation  for  the  example 
being  considered  Is  shown  In  Figure  2.8  along  with  some  addi- 
tional calculations  made  for  various  input  beam  diameters. 

For  the  conditions  described  above,  only  83  kw  of  the  initial 
400  kw  are  delivered  to  the  target.  Note,  however,  that 
by  expanding  the  Input  beam  a significant  increase  in  power 
can  be  achieved.  The  Initial  beam  size  of  10  cm  was  too 
small  to  effectively  use  the  entire  diameter  of  the  trans- 
mitting optic.  Hence,  the  beam  profile  at  the  target  was 
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Figure  2.8.  Results  of  Sample  Calculations 


larger  than  It  should  have  deen,  In  addition,  because  of 
the  mismatch  of  device  output  beam  and  beam  expander  clear 
apeivure,  a significant  power  loss  occurred.  As  the  beam 
diameter  Is  allowed  to  increase,  both  of  these  effects  are 
reduced.  Of  course,  at  the  point  where  the  beam  begins  to 
spill-over  onto  the  beam  clipper,  the  near-field  power  Is 
again  reduced  resulting  In  the  performance  fall-off  shown 
In  the  figure. 

It  Is  recognised  that  a true  performance  evaluation 
cannot  be  made  without  Including  the  propagation  effects 
between  the  transmitter  and  target  plane.  However,  the 
example  does  point  out  the  necessity  of  Including  the  effects 
of  the  optical  train  In  any  application  study,  since  It 
Indeed  has  a substantial  Impact  on  the  near-field  characteris- 
tics of  the  beam.  In  addition,  the  model  also  should  pro- 
vide some  Insight  Into  the  problems  of  Integrating  lasers 
and  pointer/tracker  systems.  It  Is  Interesting  to  note  that 
a detailed  study  of  a similar  system  using  a wave  optics 
code  resulted  In  an  additional  beam  expander  (M  » 1.3)  being 
placed  In  the  optical  train. 

The  simplified  optical  train  model  has  been  developed 
using  state-of-the-art  knowledge  of  the  contribution  of  each 
element  (windows,  mirrors,  etc.).  It  should  be  noted,  how- 
ever, that  the  model  has  not  been  verified  using  the  more 
detailed  wave  optics  codes.  Such  comparisons  would  be  use- 
ful to  not  only  verify,  but  also  to  Improve,  these  simplified 
models. 
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Sett Ion  3 
CREATION  OF  SAICQN 


One  of  the  High  Energy  liter  (HEl)  prepetition  codes 
commonly  used  by  the  User  systems  anilysts  is  the  COMBO  (or 
AT*)  code  developed  by  the  Air  Force  Weapons  laboratory 
(Reference  11).  AT*  is  the  propagation  model  documented  in 
Reference  11.  COMBO  is  a code  which  combines  ATM  with  the 
Air  Force  HEl  weight  and  volume  projections  (Reference  12). 
COMBO  also  includes  »ome  additional  features  such  ns  the 
ability  to  perform  reverse  calculations  in  which  the  user 
specifies  a desired  intensity  and  the  code  determines  the 
range  at  which  the  intensity  can  be  delivered  (given  the 
power)  or  the  power  required  (given  the  range).  The  primary 
asset  of  CCHBO  (or  AT*)  is  Its  ability  to  make  variable  al- 
titude calculations.  Including  propagation  paths  extending 
to  or  from  space. 

One  of  the  tasks  performed  under  the  present  contract 
was  devoted  to  modifying  COMBO.  The  objectives  were  to  re- 
duce execution  time.  Improve  computational  accuracy  and  to 
Incorporate  several  of  the  newer  developments  resulting 
from  the  other  tasks  of  this  contract.  A number  of  other 
minor  additions,  corrections  and  deletions  were  also  made. 

All  of  the  primary  modifications  are  discussed  in 
Section  3.1.  Collectively,  they  are  sufficient  to  warrant 
a new  code  namo  --  SAICOM.  A brief  user’s  guide  to  SAICOM 
is  provided  by  Section  3.2. 

3,1  MODIFICATIONS  OF  COMBO 

The  primary  modifications  are  dealt  with  in  individ- 
ual subsections  below.  The  basic  structure  and  organization 
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of  th«  cod*  remain  essentially  the  same  at  In  th«  origin*) 
version  (Roftrtnco  11).  A summary  of  the  wore  important 
aquations  is  given  in  Section  3.1.9. 


3.1.1  Improvements  to  Completion*!  Speed 

end  Accuracy 

A substantial  pert  of  the  COMBO  modification  effort 
wet  devoted  to  improving  the  speed  of  the  celculetions. 
Severe!  modifications  simply  involved  Improvements  in  cooing 
techniques.  These  included.  (1)  reducing  the  number  of 
cells  to  MEGAIft  (the  subroutine  which  provides  eltitude  de- 
pendent parameters  such  as  temperature*  pressure*  etc.); 

(2)  lumping  several  blooming  parameters  (viz.*  n*  dn/3T, 
and  Cp)  Into  e single  curvefit  ( GALT  IN);  (3)  using  e better 
search  algorithm  in  the  reverse  calculations;  and  (4)  reduc- 
ing the  subroutine  calls. 


Host  of  the  improvement  in  both  speed  and  accuracy 
was  obtained  by  developing  a completely  new  algorithm  for 
performing  the  double  Integration  required  to  compute  tne 
blooming  parameter,  NI.  The  expression  for  NI  is  of  the 
form 


NI  * C 


Tfp-j* / GU^di'd*1 
0 


where  C is  fixed  for  a given  set  of  conditions  (Reference 
13).  The  integration  paths  are  along  the  beam.  In  the 
original  code,  the  double  integral  was  evaluated  by  assun'- 
Ing  the  Integrals  to  be  constant  between  uniformly  spaced 
points  along  the  beam.  The  evaluation  was  successively 
repeated  with  reduced  spacing  until  the  results  for  consec- 
utive evaluations  agreed  with  10  percent.  This  procedure 
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was  found  to  be  inefficient  and  often  inaccurate.  This 
should  be  readily  apparent  from  Figure  3.1,  which  shows 
the  important  elements  of  the  integral  for  a particular 
case  of  a highly  focused  beam.  Note  that  the  value  of 
NI  is  dominated  by  the  conditions  within  the  last  few  per- 
cent of  the  range. 

The  new  algorithm  examines  the  local  behavior  of 
the  integrand  and  chooses  an  efficient  scheme  and  mesh 
spacing  dynamically.  It  selects  either  a locally  quadratic 
or  locally  exponential  repres enta t i on  of  the  integrand.  If 
neither  is  satisfactory,  it  reduces  the  mesh  spacing  and 
tries  again.  On  the  other  hand,  if  the  local  convergence 
exceeds  the  accuracy  requirements,  it  increases  the  mesh 
size  to  increase  speed.  Both  integrations  as  well  as  a third 
implicit  in  G are  performed  in  a single  loop. 

When  the  new  algorithm  was  used  for  the  severe  ex- 
ample illustrated  in  Figure  3.1,  the  calculation  required 
1.49  seconds;  the  original  COMBO  exceeded  a time  limit  of 
200  seconds  without  completing  the  calculation! 

As  mentioned  earlier,  an  additional  improvement  in 
computation  time  was  obtained  by  modifying  the  convergence 
•algorithm  used  in  the  reverse  propagation  calculation  where 
the  range  is  successively  adjusted  to  obtain  a desired  tar- 
get plane  intensity.  The  range  R'  for  one  iteration  is 
determined  by  modifying  the  previous  range,  R,  viz., 

R'  = KR^I/OI 

where  01  and  I are  the  desired  and  computed  intensities, 
respectively.  In  the  original  COMBO,  the  constant  K was 
always  1.0,  whereas  in  SAICOM,  it  varies  between  0.7  and 
1.0  to  accelerate  convergence.  In  addition,  the  convergence 
criterion  was  relaxed  from  1%  to  5%. 
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z (meters) 


Figure  3.1  Variations  in  Elements  of 
the  Blooming  Calculation 


3.1 .2  Beam  Qua! ity 

The  SAICOM  code  includes  a significant  innovation 
in  the  treatment  of  beam  quality.  It  essentially  replaces 
the  usual  one  parameter  concept  of  an  "m  x diffraction 
limited"  beam  by  a two  parameter  representation.  The  two 
parameters,  designated  and  m0,  are  used  to  specify  a 
gaussian  representation  of  the  focal  plane  intensity  profile 
so  as  to  approximate  the  actual  integrated  power  distribu- 
tion. The  two  parameter  model  is  inherently  capable  of  in- 
cluding effects  leading  both  to  beam  spreading  (m2)  and  to 
power  loss  due  to  wide  angle  scattering  (in.,), 

The  origin  and  interpretation  of  m,  and  m2  are  dis- 
cussed in  some  detail  in  Reference  1 and  in  Section  4.3  of 
this  report.  Their  use  in  characterizing  the  effects  of 
beam  truncation  and  obscuration  is  discussed  in  Reference  2. 
In  general,  they  are  calculated  in  the  new  optical  train 
model  which  is  an  integral  part  of  SAICOM  (see  Section  2). 

SAICOM  incorporates  the  two  parameter  approach  even 
though  the  user  is  only  required  to  stipulate  a single 
parameter,  m.  If  the  device  power  (PD)  and  beam  quality 
(m)  are  stipulated,  the  optical  train  portion  of  the  code 
converts  them  to  the  appropriate  aperture  or  transmitted 
power  (P.f)  and  beam  quality  (n^  , and  m2).  If  the  optical 
train  calculations  arc  bypassed  because  the  usci  specifies 
P-p  the  specified  beam  quality  is  assumed  to  be  pure 
beamspread . 

The  SAICOM  intensity  expression  for  a focused  beam 
is  of  the  form 

[1  - (1  - 0.8GS  m]  z/L)l  PTe"aZ 
iv  r 2 ( z ) 


I 


where 


r2(z)  * ~K  |l  - 2 1 2 + 4z2 
4L* 


^0.3183  m2x\ 


+ °TZ  * CTJZ]  "• 


The  term  in  the  numerator  that  decreases  as  z Increases  re- 
presents the  power  losses  that  may  be  thought  of  as  wide 
angle  scattering.  The  0.8G5  factor  is  present  because  the 

basic  propagation  equation  used  in  the  code  is  based  on  the 
2 

1/e  radiance  of  an  infinite  Gaussian  beam.  Since  and  mg 
are  used  to  account  for  the  effect  of  different  beam  pro- 
files as  well  as  phase  perturbation , these  terms  may  be  not 
the  unity  even  for  perfect  beams.  This  situation  is  clari- 
fied in  Table  3.1.  In  this  table  the  Ml  and  M2  terms  are 
the  parameters  computed  in  the  optical  train  subroutine 
while  m is  the  user  specified  beam  quality  for  the  total 
system.  Note  that  if  the  optical  train  is  bypassed  (the 
user  specifies  Py)  and  m~l,  m^  and  mg  are  0.89  and  1.29 
respectively  for  a uniform  beam  and  not  unity  because  the 
basic  propagation  equation  is  based  on  an  i nf ini te  gaussian 
model.  But  if  the  optical  train  is  not  bypassed,  mp 
Ml=0.89  and  mg-M2=1.29  because  the  basic  profile  used  in 
the  optical  train  subroutine  is  uniform. 


3.1.3  The  Optical  Train  Model 

The  development  of  an  improved  optical  train  model 
was  the  prime  objective  of  Task  5 and  is  discussed  in  de- 
tail in  Section  2 of  this  report.  The  model  has  been  in- 
corporated into  SAICOM  as  subroutine  OPTRAIN.  It  computes 
power  losses,  jitter,  and  beam  quality  degradation  associ- 
ated with  transmission  of  the  beam  from  the  laser  through 


the  final  t smltting  optic.  It  accepts  a quite  general 
speclflcati  of  the  optical  train  through  a list  of  input 
parameters.  STRAIN  accounts  for  beam  truncation  and  ob- 
scuration and  for  various  contributions  to  phase  pertur- 
bations. The  output  Includes  the  two  beam  quality  para- 
meters mj  and  mg  described  above. 

' • nmIN  can  be  bypassed  simply  by  specifying  the 
transmitted  power,  Py,  rather  than  the  device  power  In 
the  i put.  In  this  case,  beam  quality  is  based  on  wave- 
length scaling  as  usual. 

3.1.4  Multi-Line  Propagation 

Since  many  of  the  lasers  of  interest  for  high  energy 
applications  emit  several  lines  simultaneously,  it  is  de- 
sirable to  include  these  multi-line  characteristics  in  the 
simplified  propagation  codes. 

The  approach  used  in  SAICOM  is  to  compute  the  amount 
of  energy  absorbed  by  the  air  on  a line-by-line  basis.  This 
means  that  p | - Q j"R  d ( z 


' Pj  rR  ' - 

>j  irLex|)|-/  « (z'")<lz,u 
t rT  7)  1 


a(Z"  )exp 


/ Z" 

/ 


)dZ"J  in  the  integrand  of  the 


expression  for  the  blooming  parameter  becomes 


( Z" 

2J  p^:  «a1<2")exp|-y  c.e1(z»')dZ“‘j 
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where 


a#j  ■ absorption  coefficient  of  the  I**1  line 
ae^  » extinction  coefficient  of  the  1***  line. 


In  SAICOM,  the  3.8  pm  OF  propagation  has  been  extended 
to  include  three  distinct  families  of  lines.  These  encompass 
ten  lines  of  the  DF  laser,  and  represent  80%  of  the  output 
of  the  BUI  device.  With  the  basic  modification  having  now 
been  made,  it  is  straightforward  to  include  multi-line  HF 
or  CO  propagation. 

The  derivation  of  the  DF  absorption  coefficients  is 
given  in  Appendix  A and  is  based  on  the  BDL  spectral  data. 

The  resulting  S A I COM  equations  are 


where 


-A^hf-Ag+A-h) 

Aoe 

-h(6.738xlO‘5+9.GlxlO'9li) 

1 + (0.9446  RH- 1 ) e 

for  groups  1 and  2 
1 , for  group  3 

10.047,  10.082,  and  9.948  for  group  1,  2 and 
3,  respectively 

6.38RxlO"4,  G,9G2xlO'4,  and  2.199>:10'4 
2 . 25x1 0“ 8 , 1 . 55x1 0"8,  and  -1.96xl0'8 
altitude  in  meters. 


These  equations  imply  that  at  sea  level  a varies  between 
0 and  0.0409  km’1  as  relative  humidity  varies  between  0 and 
100%,  Op  between  0 and  0.0395  km’1,  and  ot^  is  a constant 
0.0478  km'1 . 
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The  above  equations  are  only  valid  to  around  h*13  km; 
for  the  higher  altitude  the  original  COMBO  coefficients  are 
used.  Unfortunately , a discontinuity  existed  at  the  13  km 
altitude;  so  transition  equations  were  created  to  connect 
the  multi-line  data  below  13  km  with  the  single-line  data 
above  that  point.  Figures  3.2,  3.3,  3.4,  and  3.5  plot  the 

V a2 * a3 * and  “resultant  e<;uat^!,s  where 

* al  * °2  * *3 
“resultant  = 3 


The  scattering  equations  are  the  same  as  those  used 
in  the  original  COMBO  (Reference  li). 


3.1.5  Modification  of  Turbulence  and  Jitter 
to  Include  Short  and  Long  Term  Effects 


During  the  code  comparison  work  in  Task  1 of  this 
contract,  it  was  noted  that  one  of  the  useful  features  of 
ESP  was  the  flexibility  of  including  beam  jitter  in  the 
blooming  calculation  if  the  frequency  were  high  enough  and 
excluding  it  otherwise  (References  1 and  14).  This  feature 
has  been  included  in  SAICOM.  The  output  jitter  factors  in 
SAICOM  are  high-frequency  jitter  (oj,^)  and  low-frequency 
jitter  (aj[_p)*  The  former  is  used  in  the  calculation  of 
the  blooming  parameter  and  both  are  combined  to  give  the 
total  beamspread  due  to  jitter.  Similarly,  the  turbulence 
calculations  are  expanded  to  include  the  total,  long-term 
turbulence  (oTLT)  and  the  high-frequency,  short-term  tur- 
bulence (ojS^.),  The  turbulence  equations  in  SAICOM  are 
based  on  work  of  Yura  (References  15  and  16)  whereas  the 
COMBO  code  was  based  on  Fried's  work  (Reference  17),  but 
the  results  differ  only  slightly  in  the  regions  of  interest. 
The  equations  used  are 

I 

jt 
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RESULTANT 


h(  km) 


Figure  3.5. 


^RESULTANT  as  4 Function  of  Altitude 


*» 


and 


"TIT 


where 


1 .439  / 

T7*  J 

LO 


c*  tt)<n 


3/S 


k - 2tt  / X 

0 * diameter  of  the  exit  . pc;  t. u . 


These  equations  arc  general i ia  1 1 on  s of  the  co-altitude  ex- 
pression developed  in  Reference  16.  In  the  blooming  ealeu 
lation  ojSj  is  used  instead  of  the  long-term  expression, 
but  it  should  he  noted  that  some  research  indicates  that 
no  turhulerce  should  be  included  in  the  hearisprrari  before 
the  blooming  calculation  (Reference  13), 

3.1.6  Far-Field  Intensity  Distribution 

In  COMBO  the  standard  calculation  of  tic  average 

9 

intensity  is  based  on  the  1/e  intensity  radius  in  the  far 
field  with  a few  other  radii  available  fur  user  spec  if i cat 
(e.g.,  1/e  and  1/2  intensity  points).  Since  SA1C0M  is  bas 

on  a gaussicn  approximation  to  the  actual  far-field  distri 

2 

bution,  it  is  a simple  matter  to  use  the  1/e  results  to 
describe  the  complete  far- field  profile.  In  addition  to 
I(l/e?),  the  peak  intensity  I ^-2. 31?.\I  ( 1/c2) , is  printed 
along  with  contour  data.  These  contour  data  assume  a 
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bivariant  normal  profile  of  the  beam  due  to  blooming  dis- 
tortion as  developed  In  Reference  ?.  For  each  102  contour 
(t*Ylpi  va0.9i  0,8,  «.,,  0.1)  the  area  ^4tMn  that  contour. 


A(I>yU  - ■ir'oSjyjd 


the  total  power  within  that  area, 


P(I>Ip)  « P(1*y) 

and  the  corresponding  average  intensity  are  printed.  In 
addition,  the  beam  displacement  due  to  non-linear  bending 
into  the  wind  and  the  eccentricity  of  the  resulting  ellptl- 
cal  proFile  are  computed  and  printed.  • — 

3.1.7  Calculation  of  Optimum  Power 

One  of  the  more  useful  bits  of  information  that 
comes  from  laser  effectiveness  studies  is  specification  of 
optimum  power  based  on  maximizing  energy  on  the  target 
(Reference  10).  Using  the  Gebhardt  and  Smith  curve-fit  for 
the  intensity  reduction  due  to  blooming  (Reference  20),  the 
far-field  intensity  is  related  to  transmitted  power,  by 

I - ,PT  tREL(N) 

where  N = bP^.  and  a and  b are  constants  (independent  of  power). 
Setting  bl/UP^O  leads  to 


Sln(IRE^)  =-3lnN 


as  the  optimum  condition.  For  the  Gebhardt  and  Smith  curve. 
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this  condition  is  satisfied  when  N=5.54.  Therefore,  the 
optimum  power  is  given  by 

5 . 5<  PT 

p s L 

optimum  N 


where  N is  the  value  of  the  blooming  parameter  corresponding 
to  the  power  P^.  The  corresponding  maximum  intensity  is 


I = 1 . 205  i 

Optimum  “ N KIP 


(KIP  is  the  computer  variable  represent 
be  noted  that  the  optimization  does  not 


n g I ^ j- . ) . It  should 
i n c 1 u <!  o the  p o s s ibis 


effects  associated  with  the  optical,  trai n , i.c.,  it  does 
not  involve  an  iterative  recall  of  OPTRAIN. 


SAICOM  prints  the  optimum  power  as  derived  above 
unless  the  corresponding  intensity  exceeds  the  air  breakdown 
threshold.  In  that  case,  it  prints  the  power  at  which  break- 
down would  occur  as  the  optimum  power.  The  breakdown  power 
level  is  estimated  on  th.  basis  of  peak  intensity  without 
blooming. 


3.1.8  Miscellaneous  Modifications 

Several  relatively  minor  additions,  deletions  and 
corrections  are  described  below. 

3 . 1 . 8 . 1 Variable  Tu r I.) u 1 o n r q in  Blooming  Computations 

In  COMBO  the  turbulence  is  not  varied  as  the  bloom- 
ing along  the  beampath  is  computed;  instead  only  the  final 
value  is  used.  In  SAICOM  the  distributive  nature  of  tur- 
bulence is  included  in  the  integral  calculation  of  the 
blooming  parameter. 
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3. 1.8. 2  Extinction  Calculations  in  Blooming  Integral 


In  COMBO  the  linear  power  losses  due  to  absorption 
and  scattering  are  computed  using  exp[-a  (h)z]  in  the  inte- 

v* 

grand  whereas  the  expression  should  have  been 

z 

exp[-  f ac(h)dz]. 

0 

This  error  causes  an  over-estimate  of  the  extinction  losses 
when  firing  downward  and  an  under-estimate  when  firing 
upward.  The  appropriate  correction  has  been  made  in  SAICOM. 

3 . 1 . 1 . 3 Power  Variation 

The  reverse  calculation  option  in  COMBO  involving 
calculation  of  the  required  power  to  deliver  a desired  in- 
tensity at  a given  range  has  been  replaced  in  SAICOM.  Now, 
whenever  the  intensity  is  calculated  in  a normal  propagation 
run,  the  power  is  automatically  varied  between  one-tenth  the 
specified  value  and  twice  the  original  power.  The  output 
then  lists  the  average  intensity  within  the  1 / a ^ intensity 
area,  the  area  over  which  the  intensity  is  above  a specified 
minimum  value,  and  the  total  power  available  within  that 
area.  The  equations  used  in  these  computations  were  dis- 
cussed in  Section  3.1.6.  The  user  is  allowed  to  stipulate 
the  minimum  intensity  of  interest  through  the  input  para- 
meter 01.  If  the  user  doe:;  not  specify  the  desired  minimum 
intensity  contour,  the  program  automatically  selects 
10  kw/cni2. 

3.1 .8.4  Deletion  of  Weight  and  Volume  Calculations 

The  technology  projection  portion  of  COMBO  is  an 
interesting  but  seldom  used  facet  of  the  code.  Therefore, 
it  has  been  eliminated. 
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3.1.9  Summary  of  Equations 

The  basic  equations  used  in  SAlCOfi  are  as  follows 

9 

Average  intensity  in  the  1/e  intensity  radius 


f. 

-E /"“ ci'z')dz' 
1-0-0. 865  m,f)l  PT  k„e  1 0 


where 


= intensity  reduction  due  tn  blooming 
(som. times  written  as  KIP) 

« f(M) 


r'(Z)  = k2rQ  + k^|  |L-/|  + /!z2  | y r-p 


/ 0.  318?!,uX  \ 2 


2 2 
+ aTLT  + aj 


aJHF  + aJLF 


r > 


L:151  ! f c2  m 

*'/*!  c"  (Z) 


The  blooming  parameter,  N,  is  computed  by 


L Z 1 

6.54  T 1 /*  |l-(l  .8r.5m.rl  F(Z,,)G(Z,,)dZ,,dZ 

~Tt  j rlT‘7  j J ^-2 

J y u / 7 ••  i f 7 " i 


where 


F'Z">  ■ PtZ)  ^ “Al'2")6  °' 


ZH 


* * r*  0 9 /0.318m.,A' 

r2(Z”)  * Vo  + k3^f  |l-Z|Z  ♦ 4ZZ  ( g ~ "t 


+ aTST  + aJHF 


aTST  “ 1 ’ 0,37  (v) 


1.J39 
“..TV?  I 
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u 

Cjj  (Z)dz 

0 
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fS(z“) 

G(z")  = iTTzlrTF(z"  )cp(z") 


3.2  A USER'S  GUIDE  TO  SAICOM 

A brief  description  of  the  key  features  of  the 
SAICOM  code  is  given  along  with  a discussion  of  how  to 
operate  the  code.  Users  familiar  with  COMBO  (Reference  11) 
should  have  little  trouble  with  SAICOM  since  many  of  the 
subroutines  are  unchanged  and  the  same  modular  format  has 
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been  retained.  The  basic  structure  of  the  executive  routine 
(SAICOM)  is  shown  in  Figure  3.6.  A substantial  portion  of 
this  routine  is  devoted  to  initializing  the  input  variables. 
As  with  COMBO,  most  of  the  input  parameters  are  defined 
through  a default  namelist  which  requires  the  user  to  specify 
only  those  parameters  which  will  differ  from  the  default 
values.  However,  it  should  be  noted  that  oven  if  the  propa- 
gation and  optical  train  namelists  are  not  used,  namelist 
cards  are  required  which  read  $ I M P U T $ and  $CASES  $,  respec- 
tively. The  optical  train  is  called  if  the  user  specified 
PD;  if  PT  is  given,  this  subroutine  is  bypassed.  A subtle 
point  worthy  of  note  is  that,  when  PD  is  given,  SIGMF  is 
assumed  to  rover  to  the  device  jitter  only,  and  tin's  term 
is  combined  with  the  four  la  jitters  (Til IP.,  TIIES,  TIIAA,  and 
T H 3 V ) given  in  n a :n  e 1 i s t C A St  to  d t ermine  l ! t d h i gh-f  requeue  y 
jitter  o':  the  total  system. 


Iterate  for  NRUNS 


Figure  3.6  SAICOH  Flow  Diagram 


Following  these  Initialization  activities,  control 
Is  transferred  to  ATM,  which  handles  the  propagation  calcu- 
lations; control  returns  to  SAICOM  only  after  the  calcula- 
tions are  completed  and  the  results  printed.  Then  SAICOM 
cycles  through  the  same  process  until  all  of  the  specified 
cases  have  been  completed. 

In  Table  3.2  a complete  list  of  the  SAICOM  subroutines 
Is  given  along  with  brief  descriptions  of  their  functions. 
Three  of  the  subroutines,  OPTRAIN,  ATM,  and  BB,  are  dis- 
cussed In  further  detail  since  they  are  the  more  important 
and  complex  portions  of  the  program. 

3.2.1  OPTRAIN 

The  calculation  procedure  within  the  Optical  Train 
Model  (OPTRAIN)  is  schematically  shown  in  Figure  3.7.  All 
inputs  to  the  subroutine  are  passed  through  a labeled  COMMON/ 
OPTIN/sta tement  and  must  be  defined  by  the  user  prior  to 
calling  the  subroutine.  Following  some  initial  calculations 
characterizing  the  input  beam,  the  effect  that  each  compo- 
nent within  the  optical  train  has  on  the  transmitted  power 
and  beam  quality  is  evaluated  in  a step  by  step  manner. 
Although  the  specific  elements  within  the  optical  train 
are  fixed,  i.e.,  aerodynamic  window,  beam  clipper,  mirrors, 
beam  expander,  and  exit  aperture,  the  user  is  allowed  to 
specify,  through  the  input  parameters,  different  types  of 
individual  components  such  as  cooled  or  uncoolcd  mirrors, 
on-axis  or  off-axis  beam  expanders,  and  open  port  or 
material  window  exit  apertures. 

The  output  from  the  model  is  returned  to  the  calling 
program  via  the  labeled  COMMO.'I/OPTOUT/  and  consists  of  (1) 
the  diameter  of  the  laser  beam  at  the  exit  aperture,  (2) 
the  power  in  the  beam,  (3)  the  total  beam  jitter  and  (4) 
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NAME 

OPTRAIN 

ATM 

AA 

Ati 

BA 

BB 

ALFSET 

ALFFAC 

ALFAAB 

ALFAEX 

CNSQ 

KI 

GAITIN 
MEGA  I R 
VX 


DESCRIPTION 

Calculation  of  Power  Loss  and  Beam 
Degradations  in  the  Optical  Train. 

Initialization  of  Propagation  Para* 
meters  and  Iteration  of  Reverse 
Cal  cul a Lions. 

Calculations  for  Laser  and  Target 
A hove  100  km. 

C a 1 c u 1 a v.  ion  f o r La  s e r A !-•  o ' ’ 100  km 
and  T r . • g e t R c 1 o v,»  100  k in . 

Calculation  for  Laser  Below  100  km  ' 
and  Target  Above  i 00  km. 

Calculations  for  Laser  and  Target 
Below  100  km. 

Sets  Power  Ratios  for  Multi-Line. 
Propayati on . 

-ot  Z -a  Z 

Calculations  of  a^e  e and  e e . 

Single  Line  Absorption  Coefficients. 

Single  Line  Scattering  Calculation. 

Atmospheric  Structure  Constant. 

Computer  Reduction  in  Intensity  Due 
to  Blooming. 

Altitude  Dependent  Parameters  Used 
in  Blooming  Integral. 

Pressure,  Density  and  Temperature  as 
a Function  of  Altitude. 

Calculation  of  Total  Crosswind. 


3-23 


Obscure 


two  parameters  (m^  and  m^}  characterizing  the  optical  quality 
of  the  beam.  The  beam  quality  parameter  Is  Internally 
adjusted  to  account  for  any  defocuslng  error  Induced  by  the 
high  power  mirrors  In  the  optical  train. 

3.2.2  ATM 

Part  of  the  Initialization  process  In  ATM  Is  to 
determine  If  the  user  has  specified  a negative  wavelength 
(e.g.,  -1.06E-5  or  -3.8E-6)  which  means  the  user  Is  speci- 
fying his  own  single  line  absorption  and  scattering 
coef  f 1 dents 

aA  » ABE  * EXP  (-ABC  * II ) 

a$  - SCL  * EXP  (-SCE  * K) 

or  his  own  pownr  ratios  for  the  multi-line  3.0  urn  propaoa- 
tlon  (see  Figure  3.8),  Note  that  single-line  3.8  pm  propa- 
gation is  available  only  if  the  user  specifies  his  own  co- 
efficients. Also  note  that  if  one  wishes  to  alter  the 
power  ratio,  FLINE(I),  I = 1,  2,  3 without  altering  the 
built-in  multi-line  coefficients,  simply  set  IAMDA  equal 
to  -3.8E-6  and  leave  the  absorption  and  scattering  coeffici- 
ents off  namelist  USPEC.  The  three  power  ratios  need  not 
total  unity  since  the  program  is  self-normalizing.  The 
heart  of  ATM  is  the  call  to  AA , AB,  BA,  or  BB  where  the 
intensities  for  the  pathlengths  of  interest  are  calculated. 

If  AB  or  BA  are  called,  they  in  turn  cal'  BB  for  those 
portions  of  the  propagation  path  that  lie  below  100  km. 
Followino  the  return  to  ATM,  the  beam  displacement,  eccen- 
tricity, contour  data,  breakdown  power,  optimum  power,  and 
effects  of  oower  variations  are  determined  for  straiqht 
propagation  unless  the  propagation  path  is  totally  exo- 
atmosoheric  (AA),  or  the  blooming  parameter  NI  is  less 
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than  10*®,  It  should  be  noted  that  when  the  power  Is  varied 
to  compute  the  Impact  of  propagating  a beam  that  Is  more  or 
less  powerful  than  the  specified  level,  the  beam  quality  Is 
not  altered;  this  elves  some  error  to  the  results  since 
bean  quality  Is  affected  by  the  total  power  transmitted 
through  the  optical  train. 

If  the  reverse  propagation  option  (ITYPE  » 5)  has 
been  chosen,  the  calculations  are  repeated  using  different 
target  ranges  ( « focal  ranges  for  ITYPE  « S)  until  the 
calculated  Intensity  Is  within  5?  of  the  desired  value. 

If  convergence  is  not  achieved  within  20  iterations,  the 
program  returns  control  to  $AIC0M  and  goes  to  the  next  case. 

3.2.3  B!i 

! 

1 In  Figure  3.9,  the  logic  flov;  for  the  PC  subroutine 

is  shorn,  all  cases  except  AA  use  this  subroutine  to  compute, 
the  atmospheric  effects  on  the  high- power  beam  propagation. 

The  subroutine  performs  calculations  at  discrete  points 

along  the  beampath,  For  calculations  ef  the  linear  propa- 
gation effects,  the  pathlergth  intervals  are  fixed  (NSTEP  - 
30  is  the  default  value)  while  in  the  blooming  loop  the 
intervals  are  increased  or  decreased  as  required  vo  minimize 
time  and  increase  accuracy  of  the  calculation.  A series  of 
linear  or  exponential  extrapolations  for  the  ith  interval 
based  on  the  i**1  - 1 interval  are  made  and  compared  with 
the  actual  calculation  for  the  ivi  interval.  If  the  compari- 
son is  unfavorable,  the  interval  is  decreased.  If  the  com- 
parison is  favorable,  the  contribution  from  the  interval  is 
incorporated  into  the  integral, if  the  comparison  is  unnecessarily 
good,  the  interval  is  increased.  This  procedure  replaces 
the  time-consuming  multiple  iterations  used  in  COMBO. 


3.2.4  Input  Data 

In  Table  3,  3 a list  of  the  input  parameters  and 
their  functions  Is  given.  This  reference  list  does  not 


Table  3.3  Input  Parameters 


REQUIRED 

ITYPE 


rn  or  n 


n 


L 

| 


If  ITYPE  » 5 reverse  calculation  is 
called  with  01  and  PD  or  PT  specific. 
All  other  entries  will  result  in 
stra  i gh  l prop  aye  i ion. 

User  may  stipulate  either  device 
power  or  tolosci  -a  ••veer ; i ? I'T  is 
specifier:  -,s  ! s,v. 

Prsirod  intensity  e : i TYPE  - F ; othcr- 
v.*  iso  it  is  too  i ' i r. ! .iiiiu  intensity  con- 
tour in  the  v,.ri.:t»u  power  calculation. 

Range  from  User  tu  Lhe  target  if  ITYPE 
i Si  if  II '.VC  = 5*  L is  a first  guess 
for  the  range. 


OPTION' 

Beam 

CHI 


DEPUG 

II 

GL 

HTM 

HOM 


MAMS  LI  SI  /IV.  HIT  / 

0 ❖ Collimated*  1 ❖ Focused. 

Horizontal  projection  of  beampath  with 
respect  to  the  direction  of  platform 
movement. 

1 ❖ no  extra  diagnostic  statement., 

2 ❖extra  diagnostic  statement. 

Focal  length. 

6round  level. 

Target  altitude. 

Laser  altitude. 


4 
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Table  3.3  (Continued) 


l 

l 


! 


O.P_U«HMt  - MMgttST;  INPUT/  (Continued) 

LAND  A Wavelength 

3 » 8x10  , 5.0x10  ® , and  1 . 0 6 x 1 0 ~ ® are 

tho  acceptable*  values.  A negative  sign 
in  front  is  a flag  that  the  user  will 
specify  the  multi-line  power  ratios  or 
the  single  line  absorption  and  scatter- 
ing coefficients. 


Seam  qual 1 ty. 


NSTEP 

Beampath  Intervals  in  the  linear  calcu- 
lation. 

OMEGA 

Slowing  rate. 

OUT 

1 *>  standard  output,  2 ->  extra  output. 

PHI 

' V r[  ^ 7 angle  with  res  poet  to  zenith. 
If  i*Hi  is  specified  Him  should  not  be 
0 i v e n a n ti  v i e e versa. 

PROP 

1 1 / e ^ 

? •; 1 '} '} !- . *•  1 ve  , t r u n e a t 0 d G a 0 s s i a n 

3 ,*  i u f i n i t e G au s s i a n . 

Rlt 

a e ; a t i v 0 humidity. 

Ri 

Radius  of  obscuration. 

RO 

Outer  radius  of  exit  aperture. 

SIGMF 

Hi  gh-f  requency  turbulence. 

VP 

P 1 a t.  r 0 mi  veloci  ty. 

VXP 

Total  crosswind  across  beam.  If  it  is 
specified  VP,  OMfGA  and  natural  wind 
are  inoperative. 

WTHR 

-1  ^high  turbulence,  low  natural  wind 

0 f*  nominal  turbulence  and  natural  wind 

1 ^ low  turbulence,  high  natural  wind. 
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Table  3.3  (Continued) 


OPT  I ON  Al.  - 

DB 

CPS 

PL 

XRI 

DIOVI 

DLI 

on L TAN 
RREF 

Xi..C 

DC 

N 

REFL 

SFAB 

ITYPE' 

XLMIR 

XMAG 

ITYPEB 


XL  FA 
BETA 
ASAB 
ITYPEA 


NAMELI  ST/  CASE/  ..  . 

Beam  diameter  coming  from  the  laser. 
Obscuration  ratio. 

Pulse  length. 

Phase  front  curvature. 

Peak-to-peak  intensity  fluctuation. 
Scale  s i e o f flue  t u a Lion  \i  r t b e am  size 
Acrov.’i  ndow . 

Aero  w i n r!  o >•/  c!  o n s i 1:  y c!  i v i c!  <•  b y a m b lent 

dens i ty . 

Distance:  from  laser  exit  to  beam  clippe 
Clipper  diameter. 

Number  of  mirrors  in  optical  train. 

Reflectivity  of  each  mirror. 

Mirror  fabrication  error. 

Mirror  type:  1 & cooled,  2 * uncooled. 

Distance  between  mirrors. 

Telescope  m a g n i f i c a t i o n . 

Telescope  type:  1 =?•  on-axis, 

2 off-axis. 

Exit  aperture  width. 

Exit  aperture  absorption  coefficient. 

Strut  area /beam  a tea. 

Type  of  exit,  aperture:  1 ->  material 

window,  2 ->  open  port. 
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Include  the  first  two  cards  which  must  precede  each  set  of 
runs.  The  first  uses  12  format  in  columns  1 and  2 to  state 
the  number  of  cases  (NRUN)  which  will  be  run,  the  second 
uses  free-field  format  to  allow  the  user  to  title  each  of 
the  runs.  The  first  card  appears  only  once  while  the 
title  card  must  precede  each  new  case.  Hence,  the  data 


stream  is  as  follows: 

NRUN  CARD 
Title 
Requi red 

Inputs 

i 

Case  1 

Namel i st 

INPUT 

Namel i st 

CASE 

Namel i st 

USPEC  (if 

r e q u i r e d ) 

Title 

Required 

Inputs 

) 

Name! i st 

INPUT 

/ 

Case  2 

Namel  i s t 

CASE 

f 

Hamel i st 

USPEC  (if 

required)  j 

Table  3, A 

provides 

a convenient 

work  sheet  format 

that  may  be  used 

to  specify 

each  case. 

It  includes  a list 

of  the  default  values  and  units.  Those  familiar  with  C0MB0 
will  recognize  a change  in  some  of  the  default  values;  the 
default  statements  now  reflect  sea  level  applications. 

Note  that  either  PHI  or  HTM  may  be  specified,  but  that  PHI 
should  always  be  given  as  a positive  angle  (measured  from 
the  zenith).  A nenative  PHI  is  a flan  that  HTM  was  specified. 

For  additional  convenience,  Table  3.5  lists  several 
of  the  important  proqram  flags  which  may  be  of  interest  to 
the  user.  A complete  program  listino  and  sample  run  is  supplied 
along  with  a sample  problem  output  in  Appendix  B. 
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Table  3.5  Flags 


FORM 

Inputs 

Negative 

Neqati ve 

Internal 

CNVRGD 

ENDED 

ERR 

FINAL 

FRACT 

LAM 

OK 


MEANING 


LAMQA  The  users  will  specify  the  multi-line 

power  ratios  or  the  single-line  absorp- 
tion and  scattering  coefficients. 

PHI  HTM  is  specified  (Default).  ■ 


TRUE  = > the  last  blooming  interval 
passed  over  the  end  point. 

TRUE  =>  this  is  the  last  blooming  interval. 

TRUE  ~>  the  reverse  calculations  have  not 
converged  in  20  iterations. 

= 1,  the  reverse  calculations  have  not 
converged  yet; 

- 2,  last  pass  throng!)  the  propagation 
c a 1 c u 1 a t i o n . ' 

Tells  the  fraction  of  the  last  blooming 
interval  that  is  within  L. 

- 1 , 10.6  ym 

= 2 , 3.8  pin  ( mul  t i-1  i no) 

= 3,  5 uni  (single-line  at  present) 

= 4,  user  will  specify  single-line 
absorption  and  scattering  coefficients. 

TRUE  =>  save  previous  values  for  multi- 
line  absorption  and  scattering  calcula- 
tions. 
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Table  3.S  (Continued) 


'I 


l 

i 

| 


{ 

|. 

t 

i 

i 


i 


FORM 

MEANING 

Internal 

(Continued) 

TYPE 

* 1,  AA 

* 2,  BB 

- 3,  AB 

» 4,  BA. 

UNIFRM  TRUE  ~>  last  two  blooming  intervals 

were  the  same. 


j. 

i 

E.  j 


•-H 

■i 

1 


J 


i 
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Section  4 

SUMMARY  OF  TASKS  I - 1 V 


The  purpose  of  this  section  is  to  briefly  review 
the  objectives  and  results  of  the  original  four  tasks. 
Specific  details  may  be  found  in  the  first  and  second 
interim  reports  (References  1 and  2).  The  primary 
objective  here  is  to  highlight  some  of  the  issues  which 
arose  during  the  course  of  the  work  and  to  indicate  some 
of  their  implications  as  they  relate  to  propagation 
model ing  in  general . 


4.1  REVIEW  OF  TASK  I.  CODE  COMPARISONS 

The  objective  of  Task  I was  to  quantitatively 
compare  t'ie  results  obtained  from  three  propagation 
codes  over  a parameter  space  representative  of  realistic 
engagement  conditions.  The  codes  were:  NOLEC  (Naval 

Surface  Weapons  Center);  COMBO  (Air  Force  Weapons 
Laboratory);  ESP-I  (United  Aircraft  Research  Laboratories). 
The  comparisons  were  presented  at  the  First  DOD  Conference 
on  High  Energy  Laser  Technology  on  3 October  1974  in 
San  Diego  (Reference  14).  An  unclassified  version  of  this 
paper  is  included  as  an  appendix  in  the  first  interim 
report  (Reference  1).  The  major  conclusions  from  these 
comparisons  were  that  all  three  codes  give  comparable 
results  for  situations  in  which  thermal  blooming  is  not 
severe,  but  that  ESP-I  gives  intensities  up  to  an  order  of 
magnitude  or  more  higher  than  the  other  two  codes  when 
blooming  is  significant.  NOLEC  and  COMBO  agree  reasonably 
well  under  most  conditions. 
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A secondary  objective  of  Task  I was  to  ascertain 
the  reasons  for  major  discrepancies  between  the  codes. 

The  effort  was  concentrated  on  COMflO  and  ESP-I  since  there 
Is  considerable  dispailty  in  their  results  even  though 
they  both  use  essentially  the  same  approach  for  calculating 
thermal  blooming.  NOLEC  uses  an  entirely  different 
approach.  In  reviewing  the  theoretical  basis  of  the  codes, 
several  Issues  were  raised  whose  significance  goes  beyond 
the  immediate  objectives  of  the  present  work. 

The  blooming  subroutines  in  both  COMBO  and  ESP-I 
are  based  on  the  theoretical  and  experimental  work  of 
Gebhardt  and  Smith  (Reference  20).  The  essence  of  the 
basic  approach  is  the  representation  of  the  ratio  of 
the  bloomed  and  unbloomed  focal  plane  peak  intensities, 
*REL’  as  a funct',on  a single  "blooming  parameter,"  N. 

N is  a theoretically  derived,  dimensionless,  similarity 
parameter.  However,  the  relationship  between  IREL  and  N 
was  obtained  experimentally.  Both  cOMBO  and  ESP-I  use 
the  same  empirical  relationship  for  IREE  (N).  They  yield 
different  results  in  blooming  calculations  becaueo  they 
compute  N differently.  The  differences  in  the  . ’’ipes 
for  calculating  N were  identified  explicity  in  thr  First 
Interim  Report  (Reference  1).  The  question  as  to  v/hich 
method  is  more  "correct"  could  not  be  resolved  on 
theoretical  grounds  for  several  reasons. 

1.  An  approximation  to  the  theoretically  derived 
expression  for  N was  used  in  correlating  the  experimental 
data.  The  approximation  was  not  generally  valid  even  for 
the  conditions  of  the  experiments.  It  is  even  more 
questionable  for  realistic  engagement-scale  conditions 
such  as  those  for  which  the  simplified  propogation  codes 
are  used. 


f j 2.  COMBO  uses  the  theoretical  expression  for  N 

l | except  for  a relatively  minor  ad,  hoc.  correction.  In 

|!  principle*  this  has  some  Important  advantages,  such  as 

I | an  Inherent  capability  to  treat  non-coaltitude  engagements. 

■ f 

[ | Unfortunately,  it  Is  not  consistent  with  the  empirical 

relationship  used  to  correlate  the  experimental  data. 

5 ■ 1 ' 

3.  ESP- I uses  the  same  expression  for  N as 
was  used  for  the  empirical  correlation.  In  that  sense, 

It  Is  at  least  consistent.  However,  it  does  not  engender 
much  confidence  In  the  extrapolation  from  small  scale 
f.\  laboratory  experiments  to  realistic  engagement  conditions. 


To  summarize  the  situation,  both  COMBO  and  ESP-I  j 

use  a semi -empiri cal  blooming  model  which  is  flawed  i 

by  a basic  inconsistency  between  its  theoretical  and 
empirical  components.  The  First  Interim  Report  implies  a 
preference  for  COMBO  simply  because  it  was  in  closer  * 

agreement  with  NOLEC.  It  should  be  emphasized,  however, 
that  this  may  be,  at  least  to  some  extent,  fortuitous. 

Furthermore,  there  does  not  exist  sufficient  data  to 
provide  convincing  validation  of  any  of  the  codes. 


I Evaluation  of  the  codes  on  theoretical  grounds 

| In  any  kind  of  absolute  sense  was  beyond  the  scope  of  the 

r 

I study,  but  two  observations  are  in  order.  First,  +he 

I consistency  problem  with  the  Gebhardt  and  Smith  blooming 

I model  does  not  necessarily  reflect  adversely  on  its 

^ basic  approach;  the  deficiencies  are  primarily  associated 

\ • with  its  implementation.  Secondly,  there,  are  no  clear 

and  obvious  theoretical  reasons  why  either  of  the  basic 
approaches  should  be  superior  in  accuracy  to  the  other. 

All  of  the  codes  incorporate  numerous  simplifying 
t : assumptions  and  approximations,  both  explicit  and  implied, 

| ; and  rely  to  some  extent  on  heuristic  arguments. 

i \ 
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4.2  REVIEW  OF  TASK  II.  BEAM  SHARE  AMO  DISPLACEMENT 

The  objective  of  this  task  Mas  to  develop  empirical 
models  for  predicting  beam  shape  distortion  and  displacement 
In  the  focal  plane  due  to  thermal  tensing.  Results  from  the 
NRl  non-linear  propagation  code  were  used  ar  he  data  base 
for  the  analysis . 

Since  beam  distortion  and  displacement  both  are 
closely  related  to  thermal  blooming,  IRa  (the  reduction 
in  peak  intensity  due  to  blooming)  was  used  as  the  primary 
correlation  parameter.  This  automatically  accounts  for 
much  of  the  dependence  on  the  independent  parameters  and 
has  the  further  advantage  of  being  readily  computed  from 
the  NRL  scaling  law. 

The  beam  displacement  model  is  represented  by 


*0,519  + 


V ( 0 . 5 1 9 ) 2 + 1.500  In  (l/IRa) 
050 


where  d is  the  displacement  of  the  peak  intensity  into  the 
wind,  af  is  the  e’^  beam  radius  (without  blooming)  in  the 
focal  plane,  a/  Is  the  e* 1 beam  radius  in  the  source  plane, 
and  n 3 af/  a^. 

The  beam  shape  model  represents  the  constant  inten- 
sity contours  as  ellipses  and  the  intensity  distribution 
as 


> 

• 

n a 

I(x,y)  a I PEAK  exp  ' 

-SO 

W- 

i * (6?; 

w 

- 
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where 


0 * 0,3  + 0.7  I 


REL 


$ » 1 


REL 


PEAK 


P .-a"  I 


REL 


This  model  automatically  gives  the  correct  peak  intensity 
•nd  conserves  total  power.  It  also  provides  an  excellent 
representation  of  the  data  in  terms  of  the  area  within  a 
given  constant  intensity  contour  and  the  eccentricity  of  the 
constant  intensity  contours. 


These  models,  of  course,  are  subject  to  the  same 
limitations  as  the  data  from  which  they  were  derived.  For 
example,  they  do  not  apply  (unless  suitably  modified)  in 
the  presence  of  jitter  or  turbulence  or  when  slewing 
is  not  in  the  same  plane  as  the  wind.  Actually,  the 
most  significant  result  of  this  task  is  that  both  distortion 
and  deflection  are  probably  too  minor  to  have  a significant 
Impact  on  system  effectiveness. 


4.3  REVIEW  OF  TASK  III.  BEAM  QUALITY 

The  objective  of  this  task  was  to  develop  an 
approach  to  modeling  the  effect  of  phase  and  amplitude 
aberrations  in  the  aperture  plane  on  the  far-field  beam 
profile.  In  more  general  terms,  it  was  desired  to  obtain 
a rational  approach  to  representing  beam  quality. 

The  recommended  approach  involves  approximating 
the  far-field  intensity  profile  (neglecting  blooming)  as  a 
Gaussian  distribution  of  the  form 
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The  parameters  and  m?  art  the  proposed  beam  quality 
parameters  and  art  determined  so  that  the  Integrated  power 
distribution  corresponding  to  the  above  expression  matches 
the  actual  distribution  at  two  specified  points*  e.g., 
r « and  r ■ 2wf,  where  w^  1$  the  beam  waist. 


This  two  parameter  representation  of  beam  quality 
Is  a significant  advance  over  the  usual  ”m  x diffraction 
limited”  concept.  It  can  account  both  for  power  loss  due 
to  wide  angle  scattering  and  for  beam  spreading.  It  should 
be  fully  adequate  to  represent  all  situations  in  which  the 
far-field  beam  profile  is  dominated  by  a central  lobe. 
Furthermore,  the  general  approach  Is  sufficiently  flexible 
that  it  can  include  effects  associated  with  beam  truncation 
and  obscuration  (see  following  section)  as  well  as  with 
other  amplitude  perturbations  and  with  phase  aberrations. 
For  example,  the  simplified  optical  train  model  (OPTRAIN) 
described  in  Section  2 provides  a method  of  estimating 
and  m.,  for  a wide  range  of  system  parameters. 


4.4  REVIEW  OF  TASK  IV.  BEAM  TRUNCATION  AND  OBSCURATION 

The  effect  of  truncation  and  obscuration  were  inves- 
tigated both  analytically  and  within  the  empirical  frame- 
work of  the  beam  quality  model  developed  in  Task  III, 
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I.t.,  the  results  were  presented  In  terms  of  the  parameters 
*1  ,n*  "2*  Oetailud  results  of  the  analysis  are  presented 
in  Reference  2.  Since  the  degree  of  truncation  and 
obscuration  Is  determined  by  the  optical  train,  the  model 

Itself  has  been  incorporated  into  the  optical  train  model 
(Section  2). 


Section  5 


CONCLUSIONS 


In  any  high  technology  program,  like  the  DOD  laser 
program,  it  is  important  that  technologists  and  analysts 
work  together  to  resolve  the  issues  and  support  the  efforts 
of  the  other  group.  During  the  early  phases  of  a new 
program  the  technologists  severly  bound  the  problem 
under  study  to  facilitate  modeling  of  complicated  processes. 
On  the  ether  hand,  the  analysts  use  extremely  simplistic 
models  for  the  elements  of  the  system  under  study  because 
they  must  address  every  aspect  of  the  problem  albeit  in  a 
naive  manner.  As  time  passes  these  two  groups  move  closer 
together  with  the  theoreticians  studying  ever  increasing 
segments  of  the  problem  and  the  systems  analysts  improving 
their  models  and  studying  the  issues  in  greater  depth. 

The  atmospheric  propagation  coding  studied  under 
this  contract  is  really  second  generation  modeling. 

Improving  the  speed  and  accuracy  of  a program,  adding 
multi-line  propagation,  and  studying  beam  shape  are  exemplary 
of  the  movement  of  HEL  system  analysts  toward  increasingly 
complex  modeling.  One  of  the  other  areas  investigated 
the  optical  train  analysis  --  is  less  well  developed. 

The  detailed  tools  ^or  modeling  the  beam  propagation 
through  the  optical  train  have  been  operational  only 
recently.  Hence,  the  simplified  codings  developed  under 
this  contract  are  first  generation  models.  The  authors 
hope  that  both  of  these  modeling  tasks  will  contribute  to 
the  continuing  covergencc  of  laser  technology  and  systems 
anal y si s . 
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APPENDIX  A 


Molecular  Absorption  of  DF  Laser  Radiation 


A.  1 INTRODUCTION 

This  appendix  was  written  by  Dr.  T.  Tuer  of  the 
Ann  Arbor  SAI  office.  This  work  was  undertaken  to  provide 
better  values  for  the  molecular  absorption  coefficients 
used  in  the  thermal  blooming  code.  These  coefficients 
are  required  as  a function  of  both  altitude  and  of 
atmospheric  water  vapor  content  for  several  "typical"  DF 
laser  line  groups.  Here,  'typical"  means  a characteristic 
laser  propagation-altitude  function  (this  will  be  clarified 
in  SectionA.3).  For  convenience  in  the  thermal  blooming  code, 
the  coefficients  are  given  as  a simple  analytical  function 
of  altitude  and  water  content,  with  different  sets  of 
coefficients  (i.e.,  one  for  each  typical  group). 

A. 2 APPROACH 

Due  to  time  constraints,  it  was  decided  to  use 
McClatchey's  DF  absorption  coefficients  ( Ref erence A . 1 ) 
rather  than  re-calculating  them  with  our  own  line-by-line 
codes.  The  ten  DDL  laser  lines  that  McClatchey  considered 
were  selected  for  this  study  (see  Table  A-l).  However, 
since  McClatchey  did  not  include  the  water  continuum 
absorption  in  his  calculations,  it  was  necessary  to  add 
thi\;  component.  A simple  computer  program  DAFT  (DF 
Altitude  Fitting  Code),  was  written  to:  (1)  evaluate 

this  continuum  component  for  the  frequencies,  altitudes 
and  model  atmospheres  desired,  (2)  add  this  to  McClatchey's 
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results,  and  (3)  plot  the  total  absorption  coefficient  as 
a function  of  altitude. 


I Table  A-l.  List  of  DF-BOL  Lines 

I 


& 

E 

DF 

Laser 

Line 

Frequency 

(cm-1) 

Relative 

Power 

Considered 
by  McClatchey 
and  Here 

[ 

P3(6) 

2546.37 

10.7 

yes 

f 

l 

P3(7) 

2570.51 

10.2 

yes 

E 

f 

P2(10) 

2580.16 

6.2 

yes 

\ 

i 

P3(6) 

2594.23 

7.7 

yes 

r 

1 

7 

P2(9) 

2605.87 

9.7 

yes 

P3<5) 

2617.41 

1.8 

yes 

P2(3) 

2631.09 

12.9 

yes 

• 

P2(7) 

2655.97 

10.0 

yes 

I 

s- 

P^  (10) 

2665.20 

3.9 

yes 

P2(6) 

2680.28 

6.2 

yes 

l 

2691.41 

9.1 

no 

I- 

£ 

P^S) 

2717.54 

6.6 

no 

P,(7) 

2743.03 

3.5 

no 

J 

i- 

■ 

P^S) 

2767  91 

2.0- 

no 

[ 

; 

The  water  cor.tiuum 

absorption 

coefficient  is  given 

L 

by  ( Reference  A. 2) . 

0*0 
n “ cs’  *pV  CN * 

P* 

w 

' 

where 

p is  the  partial  pressure  of 

water  vapor,  P is  the 

total 

pressure  and  the  * 

indicates  the  Qensi ty-equi valent- 
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.prTOi  n't  jxorp  WP75r^'^,nf)?l,l1!^wW^T:®^rrTw^?pw' 


pressures^.  This  reference  also  gives  measured  values  of 
c?  as  a function  of  frequency  at  294°  K^,  and  recommends 

**wo  n 

value  of  c«  y 8 0.12  c„  The  absorption  coefficient  n 

rl  % W S | W 

differs  from  the  usual  coefficient  k by: 
k - nU/L 

where  U is  the  absorber  thickness  and  L is  the  path  length. 
Values  for  p*,  P*  and  U were  calculated  as  a function  of 
altitude  using  data  given  by  the  Handbook  of  Geophysics 
and  Space  Environment  (Reference  A.3). 

The  attenuation  of  most  of  the  DF  laser  lines  is 
dominated  by  atmospheric  water  absorption  (i.e.,  H^O  or 
HDO  line  absorption,  or  Ho0  continuum  absorption).  In 
order  to  account  for  different  atmospheric  water  content, 
the  calculated  absorption  coefficient  was  normalized  by: 

k*  = k/p* 

Here  p*  is  a normalized  water  density-altitude  profile 

* 

which  varies  between  a value  of  p at  sea  level,  to  unity 
as  h+w.  The  quantity  pQ  is  the  sea  level  density  of  the 
model  atmosphere  under  consideration,  normalized  to  that  of 
the  midlatitude  summer  model  atmosphere.  The  form  of  p* 
was  arbitrarily  taken  to  be: 

P*  = 1 + (p*  - l)10R/? 


1 Since  the  pressures  are  generally  less  than  one 
atmosphere,  we  used  the  usual  pressure  in  place  of 
the  densi ty-equiva i ent- pressure . 

^ Since  the  temperature  dependency  ef  this  quantity  is 
uncertain,  its  value  at  294°K  was  used  throughout. 
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where  R is  the  normalized  water  density  profile  based  on 
the  data  in  Reference  A.^  (see  Figure  A^l).  An  example  of 
the  effect  of  this  scaling  is  shown  in  Figure  A-?,  where 
the  absorption  coefficient  profiles  for  the  five  different 
model  atmospheres  are  seen  to  collapse  nicely  into  a narrow 
band , 


A-4  ?/ 


— -r 


\ ^vO 
\ > 
v \ N 

Vx  \ \ 

s ' N > 

V W\ 

\ ''.A,  N 

V >v  * \ \ 
v . Nv  s \ fiN 


V \ 


Atmospheric 

Model 

Tropical 

Midiatitnde  Summer 
Midlalitutfe  Winter 
Subarctic  Summer 
Sub  rette  Winter 


X\  V /-*  - Cjh  + c2h 

Vv  \S  c,  « 0.  Of  D73 

\\  A 1 


Or,  =“  0.  C?o'  5 


X '+ 


\\  \ 

\N\  \ 

>4'  \ 

W « 


Note:  Water  vapor  density  < 
(gm/m^)  can  be  con- 
verted to  precmiinble 
centimeters  err.  :n) 
bv  dividin'*  bv  lb. 

* V.  » 


\A  V 

* \\ 

X \°V 

\ \ \ 

A \\ 

A ' \ 

\ A \ \\ 

VA  V V 


Altitude  (kml 

Figure  A-l.  Water  Concentration  Profiles  for  Five  Model  Atmospheres 
and  Analytic  Fit. 
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A. 3 DETERMINING  TYPICAL  GROUPS 

Normalized  profiles  for  all  the  water  dominated 
Ustr  lines  were  plotted  on  a common  scale,  hf  course 
the  lines  which  are  dominated  by  other  species  (e.q., 

HjjO  line  absorption),  were  not  normalized  in  this  fashion 
but  simply  plotted.  These  plots  tend  to  collect  in  throe 
general  groups.  Those  laser  lines  which  are  dominated  by 
water  at  low  altitudes*  but  happen  to  fall  on  underlying 
absorption  lines  (Type  A),  have  a knee  in  their  profile 
(see  Figure  A-3).  Those  that  do  rot  have  underlying 
lines  (Type  8),  do  not  have  a pronounced  knee  but  continue 
to  decrease  with  altitude,  a:  least  foe  *.ne  altitude  ranee 
considered  here  (see  Figure  A-4).  Lines  which  are  not 
dominated  by  water  (Type  C),  have  no  tree  but  fall  off 
less  rapidly  with  altitude  (sec  Figure  A-5). 

The  mechanism  fur  each  type  can  be  elucidated  by 
examining  the  details  of  the  absorbing  lines  and  continuum. 
Figure  A-6  shows  components  of  the  absorption  coefficient 
of  a Type  A laser  line  at  sea  level.  Notice  the  under- 
lying CH^  absorption  line  (M2),  which  will  become  important 
at  higher  altitude  as  the  absorption  by  water  diminishes. 
Figure  A- 7 shows  a Type  B laser  line  where  there  is  no 
important  underlying  absorption  lines.  Finally,  a Type  C 
laser  line  is  shown  in  Figure  A-b,  which  is  dominated  by 
NjjO  absorption  even  at  sea  level,  and  is  affected  little 
by  water. 


A. 4 RESULTS 

A least  squares  fit  was  applied  to  each  group  of 
laser  lines  (i.e..  Types  A,  B and  C).  The  analytic  t\:m 

used  was: 

In  k*  « AQ  + Ajh  + A?h2 


rvr.rucs  cut) 


Figure  A~3.  leaser  Prenerntior ,-AVitudo  Function  Typo  A:  I>or:ii»’ntod  by 
Water  at  Low  Ai;iuKV:i  but  with  Underlying  Absorption  Lines 
which  become  Important  at  High  Altitudes. 


Figure  A--L  Laser  Propagation-Altitude  Function  Type  U:  Dominated  by 
Water  at  Low  Altitudes  but  with  Underlying  Absorption  Lines 
which  become  Important  at  High  Altitudes. 
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Figure  A-7.  Contributors  to  the  Molecular  Absorption  of  the  P?(9) 

DF  Line  (Midlatitude  Summer,  Sea  Level). 

Example  of  Type  B Absorption,  Extracted  from  Ref/  A . 4 
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Freqreucy  (cirT*) 

Figure  A-8.  Contributors  to  the  Molecular  Absorption  of  the  P„(8) 

DF  Line  (Midlatitude  Summer,  Sea  Level). 

Example  of  Type  C Absorption,  Extracted  from  Ref.  * , 4 
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where  k*  is  the  normalized  absorption  coefficient  (km"*)  and  h is  the 
altitude  (km).  The  results  of  this  fit  for  each  group  are  shown  in 
Figures  A-9  through  A-ll,  and  the  determined  coefficients  (i.e. , Aq, 
Ap  Ag)  given  in  Table  A -2. 

Table  A-2.  Coefficients  of  Least  Squares  Fit 


Type 

DF  Laser  Lines 

Aq 

A! 

A2 

A 

P2(8),  P2(7),  Pj(10) 

-3. 138 

-0.  6388 

0.  0225 

B 

P3(6),  Pa(9),  P3(5),  P,(6) 

-3. 173 

-0.  6952 

0.  0155 

C 

Ps(8),  P/7),  P2(10) 

-3.  039 

-0.  2199 

-0.  00196 

For  Type  C laser  linos  this  normalized  coefficient  k*  is  identical  to  the 
un-normalized  coefficient  k (since  water  absorption  is  insignificant): 

k = k*  (for  Type  C) 

Type  A and  3 lines  must  be  un-normalized  by: 

k»p*k*  (for  Types  A & B) 

where  P*  si  + (p*  . i)  ioR/2  and:  ' “ .■  - (( * -i)  io^"'  A ; 

R = CAh  + C2h2  \ 

with  Cj  = 0.05375  and  Cg  " 0.00835. 
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Figure  A-P.  Least  Squares  Fit  to  Type  A Laser  Lines 
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